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ABSTRACT 


This  report  (Volume  I,  Empirical  Methods)  presents  the 
methods,  equations,  and  substantiating  data  for  an  empirically 
based  computer  program  for  the  rapid  and  accurate  evaluation 
of  the  aerodynamic  characteristics  of  large  aircraft  (bombers, 
tankers,  and  transports)  from  takeoff  through  landing  and 
through  the  subsonic,  transonic,  and  supersonic  speed  regimes. 
The  program  calculates  lift,  moment,  and  drag  characteristics 
at  both  low-  and  high-lift  conditions,  including  the  effects 
of  ground  proximity  during  landing  and  takeoff.  The  Input 
requires  the  configuration  geometry  and  the  aerodynamic  condi¬ 
tions  for  which  solutions  are  desired.  The  program  includes 
the  capability  of  analyzing  both  fixed-wing  and  variable-sweep¬ 
wing  configurations  as  well  as  the  aerodynamic  characteristics 
of  the  most  recent  supercritical  wing  designs.  The  accuracy  of 
the  program  is  verified  through  comparisons  of  the  predicted 
results  with  experimental  data  for  several  configurations. 
Details  of  the  input  and  output  for  this  program  along  with 
a  FORTRAN  source  deck  listing  and  sample  problem  are  contained 
in  Volume  II.  Program  User  Guide.  Although  this  program  was 
developed  to  handle  the  bomber,  tanker,  transport  class  of  air¬ 
craft,  it  is  also  applicable  to  fighter  type  aircraft  without 
maneuver  devices. 
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1.  INTRODUCTION 


The  need  for  a  computer  program  for  performing  rapid  and 
accurate  evaluation  of  the  aerodynamic  characteristics  of  an 
aircraft  is  evident  during  evaluation  of  a  preliminary  design 
configuration  when  both  time  and  geometry  definition  are 
limited.  A  quick-response  program  is  also  needed  for  performing 
preliminary  design  trade-off  studies.  This  report  presents 
documentation  for  an  empirically  based  computer  program  that 
will  predict,  with  minimum  aircraft  geometry  input  requirements, 
the  aerodynamic  characteristics  of  large  aircraft  (bombers, 
tankers,  and  transports)  from  takeoff  and  landing  through  the 
subsonic,  transonic,  and  supersonic  speed  regimes. 

Some  of  the  methods  contained  in  the  AeroModule  program 
(Reference  1)  developed  at  General  Dynamics '  Convair  Aerospace 
Division,  Fort  Worth  Operation,  are  utilized  in  the  Large- 
Aircraft  Aerodynamic  Prediction  Program.  The  AeroModule  pro¬ 
gram  was  developed  to  provide  lift  and  drag  estimates  for  use 
in  a  computerized  aircraft  design  synthesis  program.  The 
Large-Aircraft  program  has  extended  the  AeroModule  methods  to 
include  moment  calculations  along  with  improving  the  methods 
used  to  predict  lift  and  drag. 

Input  to  the  Large  Aircraft  program  requires  the  configu¬ 
ration  geometry  and  the  aerodynamic  conditions  for  which  solu¬ 
tions  are  desired.  Aircraft  geometry  is  represented  by  a 
series  of  component  bodies  and  airfoil  surfaces.  The  geometry 
input  requirements  are  minimized  by  use  of  internally  calculated 
values  where  possible.  Component  wetted  areas  can  either  be 
calculated  internally  from  other  input  or  input  directly, 
as  desired.  The  program  includes  the  capability  of  analyzing 
both  fixed-wing  and  variable-sweep-wing  configurations  along 
with  either  conventional  or  supercritical  wing  designs. 

The  output  from  the  Large  Aircraft  program  consist  of 
tabulated  Cl,  Cp,  Cm,  and  angle-of-attack  predictions  for  a 
given  Mach  number  and  altitude  or  Reynolds  number  condition. 

In  addition,  breakdown  is  given  of  each  drag  item  such  as 
drag  rise,  wave,  friction,  base,  trim,  camber,  etc.  Also, 
other  aerodynamic  parameters  such  as  lift-curve  slope,  CLmax, 
polar-shape  factors,  wing-body  aerodynamic  center,  etc.,  are 
listed.  The  output  format  is  designed  to  provide  the  program 
user  with  a  quick  scan  of  the  significant  aerodynamic  para¬ 
meters  along  with  an  in-depth  lock. 
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The  methods,  equations,  and  substantiating  data  for  the 
Large  Aircraft  program  are  presented  in  the  following  sections. 
Details  of  the  input  and  output  of  this  program  along  with  a 
FORTRAN  source  deck  listing  and  sample  problems  are  contained 
in  Volume  II,  Program  User  Guide. 

Although  this  program  was  developed  to  handle  the  bomber,  tanker, 
transport  class  of  aircraft,,  It  Is  also  applicable  to  fighter  type  air¬ 
craft  without  maneuver  devices.  The  program  Is  referred  to  throughout 
the  tpxt  as  the  Large  Aircraft  Program. 


2.  GEOMETRY 


The  Large- Aircraft  Aerodynamic  Prediction  Program  requires 
a  minimum  of  input  data  since  most  of  the  geometric  parameters 
used  in  the  aerodynamic  methods  are  calculated  internally  within 
the  program.  Some  geometric  parameters  such  as  wetted  areas  and 
mean  geometric  chords  can  be  either  generated  internally  by  the 
program  or  accepted  as  input  data.  The  conventions  and  equations 
used  by  the  Large  Aircraft  program  to  determine  the  geometric 
parameters  used  in  aerodynamic  calculations  are  described  in  this 
section. 


2 . 1  Component  Geometry 

The  basic  aircraft  geometry  is  represented  by  a  series  of 
components.  The  fuselage,  canopy,  and  stores  are  represented  by 
a  series  of  bodies;  the  nacelles  are  represented  by  another  series 
of  bodies  (open nosed) ;  and  the  wing,  tail  surfaces,  pylons,  and 
ventrals  are  represented  by  a  series  of  single-panel  airfoil  sur¬ 
faces.  For  cranked  or  complex  wing  planforms  the  wing  can  also 
be  represented  by  a  series  of  interconnected  surface  panels. 
Provision  is  also  incorporated  into  the  program  for  computing 
the  geometric  parameters  for  variable-wing-sweep  configurations. 


2.1.1  Body  Geometry 

The  minimum  geometry  input  requirements  for  the  body  compo¬ 
nents  are  length,  width,  height,  nose  length,  and  boattail  length; 
in  addition,  for  open-nosed  bodies,  the  inlet  and  exit  area  must 
be  specified.  If  the  maximum  cross-sectional  area  of  the  compo¬ 
nent,  Amax»  is  not  input,  the  value  is  calculated  by 

aMAX  *  J  (width  x  height)  (2-1) 

Also,  if  wetted  area  for  the  component  is  not  input,  a  value  is 
then  calculated.  The  wetted  area  for  close-nosed  bodies  is 
determined  by 

awet  "  +  2,5/bt(1  +  ^  ^"majP ) 

+  4(/-  /n  -  AMAX  <2"2) 
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For  open-nosed  bodies,  wetted  area  is  determined  by 


■  /bT^J  \A 4  AMAX 


(2-3) 


2.1.2  Airfoil  Surface  Geometry 

The  input  required  to  define  the  planform  for  the  airfoil 
surface  components  that  represent  surfaces  other  than  the  main 
wing  are  exposed  taper  ratio,  exposed  root  chor4  and  the  leading- 
and  trailing-edge  sweep  angles. 

The  aspect  ratio  and  the  exposed  area  of  each  component 
are  calculated  by  the  equations 


AK 


(i-A,) 


s  7T+"57T  tan 

'  v  /  g  - \  it*/ 


(2-4) 


JLCi'  s 


JEXP 


[(Cr)s  (1+Ss)] 


2  ARc 


(2-5) 


If  the  component  wetted  area  and  mean  geometric  chord  are  not 
input,  they  are  calculated  by  the  equations 

1  +  4^1  1  (2-6) 


f 


r 


SWET  =  SEXP  |_2+.1843(t/c)g  +  1 . 5268( t/c)g2 
-  .8395(t/c)fi3J 


(2-7) 


The  wetted  area  is  essentially  twice  the  exposed  area  with  a 
small  factor  to  account  for  thickness. 

For  trim  calculations  the  location  of  the  horizontal  tail 
quarter-chord  point  on  the  MAC  is  calculated  from 
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(2-8) 


V  •— »  111 

XHT 

where  by?  is  the  exposed  area  span  of  the  horizontal  tail  calcu¬ 
lated  from  bHT  “  V  ARHT  *  SEX?  The  moment  arm  of  the  tail  at 
any  angle  of  attack  is  computed  from 

/^HT  "  XH  ccs  (-0.  -OC.  )  (2-9) 

where 

XH  *  |/ (ZHT  -i  ZCG)2  +  (XHT  -  XCG )* 


tan<'WHT  +  Tp  +  <xle) 


HT 


and 


Jl 


arc tan  < 


f (2HT  ~  zcg) 
<%T  "  XCg) 


XGG  and  ZGG  are  the  longitudinal  and  vertical  locations  of  the 
moment  reference  point  and  Zj|j  is  the  vertical  location  of  the 
quarter- chord  point  on  the  horizontal  tail  Mac. 

2.1.3  Wing  Geometry 

If  the  main  wing  is  defined  as  one  panel,  the  aspect  ratio, 
taper  ratio,  leading-edge  sweep,  planform  area,  and  the  location 
of  the  wing  relative  to  the  fuselage  are  required  input  to  define 
the  planform  geometry.  The  geometry  for  a  complex  wing  planform 
is  represented  as  a  series  of  panels,  as  shown  below  in  the  sketch. 


The  total  wing  is  defined  by  up  to  ten  panels.  The  leading-edge 
location  and  chord  length  of  each  panel  edge  is  specified  along 
with  the  average  section  camber  and  thickness  of  each  panel. 
Average  values  of  thickness  and  camber  are  computed  by  the  root- 
mean-  square  equations 


t/c 


i 


(2-10) 


XN(cia>t 

SEXP 


Si 


(2-11) 


where  Sgxp  ~  Si  is  the  sum  of  the  exposed  areas  of  all  the 
panels.  Certain  aerodynamic  calculations,  such  as  wing  wave 
drag  and  lift- curve  slope,  require  the  use  of  an  "equivalent" 
trapezoidal  wing  that  approximates  the  planform  of  the  arbitrary 
wing.  The  sweep  angle  of  the  equivalent  wing  is  obtained  by 
area-weighting  the  sweep  according  to  the  equations 

N 

2  (tan/tLE>iSi 

(tan\E)e  «  ^ - —  (2-12) 

SEXP 

N 

2  (cosv\c/4>isi 

(cosAc/4)e  -  — -  (2-13) 

sExr 


N 

2  (c°sA:/2)isi 

(cosXc/2)  “  — — -  (2-14) 

'  e  c 

dEXP 


(tan^E)e 


( tan/q^g) 


SEXP 


(2-15) 
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If  the  wing  panel  wetted  areas  and  mean  geometric  chord 
are  not  input,  they  are  calculated  by  use  of  equations  similar 
to  Equations  2-6  and  2-7. 

For  a  multiple-panel  wing,  if  the  planform  area  is  not 
input,  a  value  is  calculated  by  summing  the  panel  exposed  areas 
and  adding  the  area  obtained  by  extending  the  innermost  panel 
to  the  centerline  of  the  aircraft.  In  the  case  of  a  wing  whose 
innermost  panel  represents  a  strake  with  a  large  leading-edge 
sweep  angle,  extending  this  panel  to  the  centerline  of  the  air¬ 
craft  would  result  in  an  extremely  large  planform  area.  In 
this  case  the  value  of  the  theoretical  planform  area  of  the  wing, 
ignoring  the  strake,  should  be  input.  The  aspect  ratio  is  de¬ 
fined  as 

AR  -  b2/SpLAN  (2-16) 

Lift  and  drag  parameters  are  calculated  by  use  of  the  aspect 
ratio  defined  with  the  wing  planform  area,  and  the  final  results 
are  then  referenced  to  the  reference  area,  which  is  input. 

In  most  typical  cases,  SprAM  equals 


2 . 2  Variable-Sweep  Configuration 

The  planform  for  a  variable- sweep  configuration  Is  defined 
by  a  trapesoidal  movable  panel  and  an  optional  glove  panel,  as 
shown  in  the  sketch  below.* 
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The  procedure  first  defines  the  coordinates  of  Points  1,  2, 
3,  and  4  from  the  input.  (The  input  planform  area  is  equal  to 
twice  the  area  enclosed  by  these  four  points.)  When  the  movable 
panel  is  rotated  about  the  pivot  point,  the  resulting  geometry  is 
as  sketched  below.  The  coordinates  of  Points  2,  3,  5,  and  6  are 
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then  determined.  The  planform  area  is  calculated  as  twice  the 
area  enclosed  by  these  four  points.  Since  the  tip  chord  is  assumed 
streamwise  at  the  forward  reference  sweep,  the  distance  from  Point 
2  to  the  centerline  is  the  semi- span,  b/2.  The  aspect  ratio  is 
thus  defined  as 

AR  “  / Stit  (2-16) 

The  taper  ratio  is  calculated  as 

(2-17) 


where  Cy  and  CR  are  as  defined  in  the  following  sketch. 
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The  mean  aerodynamic  chord  (MAC)  and  wetted  area  of  the 
outboard  panel  are  calculated  by  use  of  the  tip  chord,  Op,  and 
the  chord,  CRX2-  The  (MAC)calculat-ed  is  expressed  for  the  out¬ 
board  panel  as 

J  +  1  (2-18) 

1+A  J 

where  $  -  CT/CRx2.  The  wetted  area  (Awet)caicuiated  is  com¬ 

puted  using  twice  the  exposed  area  of  the  panel  and  the  thick¬ 
ness  correction  to  wetted  area  expressed  by  Equation  2-7.  These 
calculated  values  for  MAC  and  A*,^  are  compared  with  the  optional 
input  values  at  the  forward  reference  sweep  and  the  aft  reference 

o  Moon  nnn/1  ■(  t1!  nno  T  f  aoI  anl  4  n  im  i  f>  tmltmo  rl  4  on 

“  "  wwff  •*«  «-  -.V#  **»-»  •  4.  *■  vuv.  Vt4  4.VV»4.fc*  V  V-Vi  U  HV»  VQliUVO  U  1.1.1.  Wli  j  Cl  a 

might  occur  for  non- trapezoidal  planforms,  the  input  values  are 
used  and  the  incremental  differences  between  the  two  are  used  for 
interpolation  purposes  in  the  calculation  of  MAC  and  at 

intermediate  sveep  angles.  The  equations  are  as  follows: 

MAC  -  HACcalc  +<4MACi  +  (AMAC2-AMACp)  -  (2-19) 

^2'Ai 

^wet  “  ^wetca^c  +  ^^wetp  +  (^wet2"1^wet^)^^_^  (2-20) 
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where 


^MAC  -  MACcalc.  -  MAClnput 
■^Awet  *  Aweccaic#  *  A«etinput 


and  where  the  subscript  1  refers  to  the  forward  reference  sweep 
position,  and  the  subscript  2  refers  to  the  aft  reference  sweep 
position. 


The  maximum -thickness  sweep  angle,  A(t/c)max>  usec*  *n  fric¬ 
tion  drag  calculations,  is  calculated  from  the  quarter-chord 
sweep  of  the  panel  and  the  inputjV( t/c)tnax  at  f°rward  and 
aft  reference  sweep  positions.  The  equation  is 


~Ac/4  +4-yli  +  (442-4A.1)j|— jj- 

(2-21) 

where 

(Ac/ 4)calc.  -  <Ac/c)max)lnput 

The  streamwise  camber  and  thickness  of  the 
at  a  given  sweep  are  calculated  by 

outboard  panel 

^L^calc.  ^L^ref  (c'/cjref 

(2-22) 

(t^c)calc.  (t/c)ref  (C'/C) _ r 

(2-23) 

where 

r./c  _  0  5  f  cos^LE  cosAte  7 

C  / c  -  0.5  /  j  _  ; — rrv  +  -'tt  l \  > 

(2-24) 

Equation  2-24  is  the  relationship  between  the  chord  perpendicular 
to  the  mid-chord  sweep,  C’ ,  and  the  streanrwise  chord,  C.  For  a 
variable-sweep  wing,  C'  remains  constant  so  that  the  camber  and 
thickness  perpendicular  to  the  mid-panel  sweep  also  remain 
constant.  Finally,  the  outboard  panel  thickness  is  compared  with 
the  aft  reference  sweep  input  value  and,  if  the  calculated  and 
input  values  differ,  the  input  value  is  used  for  interpolation 
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purposes  in  the  calculation  of  t/c  at  intermediate  sweep  angles. 
The  equation  is  as  follows: 

I 

(t/c)  -  <t/c)calc  +  0t/c)  (2-25) 

where  (At/c)  -  (t/c)input  -  (t/c)calc  at  aft  sweep. 

The  variation  of  wing  twist  with  sweep  can  be  calculated 

from 

T -  arctan  (—^)  (2-26) 

L  tip 

where  AZ  is  the  vertical  position  of  the  leading  edge  (assuming 
the  wing  is  twisted  about  the  trailing  edge)  and  C  is  the 
streamwise  chord  at  the  tip  in  the  swept  position.  The  tip 
displacement  is  calculated  at  the  forward  reference  sweep  posi¬ 
tion  through  the  equation 

(AZ)fcip  *  (CR^)  tan  X. 

The  streamwise  chord  at  the  tip  is  calculated  from 
ctip  “  CR  “  <b/2)(tanALE  -  tan/TE) 

The  tip  displacement  is  assumed  to  be  independent  of  sweep. 

The  variation  of  wing  incidence  with  sweep  is  calculated 

from 

i  ^  iref  "  tanA/LE  *  tanAxEl^  cosAAle  (2-27) 

'  *v 

where  "WEE  -ALE  .  In  calculating  the  variations  of  twist 

and  incidence  with  wing^sweep,  it  is  assumed  that  the  wing  pivot 
is  perpendicular  to  the  wing  chord  plane. 


2 . 3  Airfoil  Section  Geometry 

Several  airfoil  section  parameters  are  used  in  the  aerodyna¬ 
mic  predictions.  These  parameters  are  generated  internally  in  the 
program  for  the  NACA  6-series  and  4-digit  airfoil  sections  along 
with  biconvex  and  supercritical  airfoil  sections.  The  procedure 
determines  the  leading-edge  radius  as  a  function  of  thickness 
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ratio,  t/c,  for  these  airfoils,  as  shown  in  Figure  1.  Hie 
distance  of  the  position  of  maximum  thickness  from  the  leading 
edge,  Xt/cmax*  is  listed  in  Table  I.  A  leading-edge  sharpness' 
parameter, Ay»  expressed  as 

Ay  -  A(t/c)  (2-28) 

is  defined  for  uncambered  airfoils,  where  A  is  a  function  of  the 
airfoil  leading-edge  geometry  (shown  plotted  in  Figure  2.2. 1-8 
of  Reference  2).  The  trailing-edge  angle  of  the  upper  surface 
of  the  airfoil  is  computed  from 

0TE  -  B(t/c)  +  C(C/(J)  (2-29) 

A,  R,  C  values  used  in  the  Large  Aircraft  program  are  listed  in 
Table  I. 

If  the  airfoil  section  cannot  be  approximated  by  one  of  the 
sections  contained  within  the  Large  Aircraft  program,  the  user 
can  Input  geometry  to  define  any  arbitrary  airfoil  section. 

Two  examples  of  the  designation  for  a  six-series  airfoil  are 
given  by? 

64-  210  and 
64A210 

The  6  for  the  first  digit  indicates  a  6-series  airfoil.  The 
second  digit  (4)  designates  the  chordwise  location  (in  tenths) 
of  the  minimum  pressure  for  the  basic  symmetric  airfoil  at  zero 
lift.  The  third  digit  (2)  designates  the  camber  design  lift 
coefficient  (in  tenths).  The  last  two  digits  (10)  designate 
the  airfoil  thickness  (in  percent).  The  letter  A  appearing  in 
some  6-digit  series  designations  indicates  that  a  modified 
thickness  and  camber  distribution  is  used. 

An  example  of  the  designation  for  a  4-digit  airfoil  is  given 

by: 


0012-34 

where  the  12  designates  the  thickness  (in  percent  chord),  the  4 
designates  the  position  of  maximum  thickness  (in  tenths),  and 
the  3  designates  the  leading-edge  radius  (3  designates  1/4  normal, 
6  designates  normal,  and  9  designates  9X  normal  leading- edge 
radius). 
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Table  1 


TABULATED  AIRFOIL  SECTION  PARAMETERS 


Airfoil  Type 

^(t/c)max 

A 

B 

C 

63-series 

.35 

22.0 

34.6 

14.8 

64-series 

.375 

21.7 

38.4 

14.8 

65 -series 

.41 

19.2 

46.4 

14.8 

66-series 

.44 

18.35 

60.2 

14.8 

63A 

.37 

22.0 

57.5 

14.05 

64A 

.39 

21.2 

59.5 

14.05 

65A 

.42 

19.2 

66.5 

14.05 

Supercritical 

.3471 

27.0 

30 

40.0 

Biconvex 

.50 

11.75 

95.0 

0.0 

00XX-62 

.2 

24.0 

50.0 

13.8 

-63 

.3 

24.0 

63.0 

13.8 

-64 

.4 

22.0 

82.8 

13.8 

-65 

.5 

20.0 

113.0 

13.8 

-66 

.6 

20.0 

153.0 

13.8 

-33 

.3 

19.0 

63.0 

13.8 

-34 

.4 

17.0 

82.8 

13.8 

-35 

.5 

15.0 

113.0 

13.8 

-93 

.3 

29.0 

63.0 

13.8 

-94 

.4 

27.0 

82.8 

13.8 

-95 

.5 

25.0 

113.0 

13.8 

3.  MINIMUM  DRAG 


The  drag  of  an  aircraft  can  be  represented  as  the  sum  of 
minimum  drag,  plus  drag  due  to  lift,  plus  drag  due  to  trim.  The 
drag  bookkeeping  system  used  in  the  Large  Aircraft  program  has 
minimum  drag  comprised  of  the  drag  items  that  are  "assumed"  to 
be  independent  of  lift,  such  as  friction,  form,  interference, 
wave,  base,  camber,  roughness  and  proturberance.  Drag  due  to 
lift  is  comprised  of  the  drag  items  that  vary  with  lift,  such 
as  induced  drag,  profile  drag  increment  due  to  lift,  and  flow 
separation  drag.  Transonic  drag  rise,  which  varies  with  lift, 
is  separated  for  bookkeeping  purposes  into  an  increment  added 
to  minimum  drag  and  an  increment  added  to  drag  due  to  lift.  In 
cases  where  the  fuselage  has  an  upswept  aft  end,  the  Increment 
in  fuselage  drag  between  an  upswept  fuselage  and  a  symmetrical 
fuselage  is  tabulated  in  the  program  output  as  a  function  of 
lift.  The  drag  buildup  does  not  include  Incremental  drag  con¬ 
tributions  due  to  propulsion  installation  such  as  spillage  drag, 
bleed,  nozzle  effects,  etc.  In  many  thrust-drag  accounting  sys¬ 
tems  the  propulsion  related  drag  increments  are  included  in  the 
propulsion  force  buildup  since  these  drag  increments  vary  with 
power  setting.  If  a  horizontal  tail  is  present  on  the  configu¬ 
ration,  the  untrimmed  lift  and  drag  is  computed  for  a  zero  tail 
deflection  condition.  The  effect  of  horizontal  tail  deflection 
for  trim  is  determined  by  computing  the  lift  and  drag  increment 
relative  to  the  zero  tail  setting. 

The  methods  used  to  determine  each  of  the  minimum  drag 
contributions  and  the  fuselage  aft-end  upsweep  drag  are  de¬ 
scribed  in  the  following  subsections.  Drag  rise,  drag  due  to 
lift,  and  trim  drag  are  discussed  in  Sections  5,  4,  and  7, 
respectively. 


3.1  Friction,  Form,  and  Interference  Drag 


A  large  part  of  the  subsonic  minimum  drag  is  comprised  of 
the  sum  of  friction,  form,  and  interference  drag  of  all  the 
aircraft  components.  The  drag  of  each  component  is  computed  as 


CD 


(Cf  * 


Awet 

SREF 


)  •  FF  •  IF 


(3-1) 


where  Cf  is  the  compressible  flat-plate  skin- friction  coeffi¬ 
cient,  A^et  *-s  the  component  wetted  area,  and  FF  and  IF  are  the 
component  form  and  interference  factors. 
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3.1.1  Friction  Drag 


The  flat-plate,  compressible,  turbulent,  skin- friction  coef¬ 
ficient  is  determined  from  the  general  equation  given  in  Refer¬ 
ence  3, 

Cf  -  JY  Cfi(RNL  •  F2>  0-2) 


where  Fi  and  F2  are  functions  of  the  frees tream  Mach  number  and 
wall  temperature.  The  incompressible  skin- friction  coefficient, 
Cf,,  is  evaluated  at  the  equivalent  Reynolds  number,  R^  F2* 
White  and  Christoph  (Reference  3)  developed  expressions  for  the 
transformation  functions  F^  and  F2  along  with  a  more  accurate 
explicit  equation,  based  on  Prandtl/Schlichting  type  relations, 
for  computing  the  incompressible,  turbulent,  flat-plate  friction 
coefficient  (Cf^)  .vith  the  following  results: 


F2  -  t1+n  f 


0,430 


H  (lo*10  rnl)2  ' 


56 


For  an  adiabatic  wall  condition. 


t  and  f  are  given  by 


t  -  Tw/Taw  “  [l  +  r 
f  -  1  +  0.044  r  M„2  t 


Using  a  recovery  factor  r  •  0.89  and  a  viscosity  power-law  expo¬ 
nent  n  ■  0.67,  recommended  in  Reference  3,  results  in  the  follow¬ 
ing  expression  for  Cf: 


Cf  -  t  f" 


0.430 


l°SlO<%L 


•  f)) 


275S 


(3-3) 


where 

t  -  [l  +  0.178  Mj2]’1 
f  -  l  +  0.03916  M  2  ■  t 

OO 

The  Reynolds  number,  Rjj^,  is  based  either  on  component  length  or 
an  admissible  surface  roughness,  whichever  produces  a  smeller 
value  of  Reynolds  number,  as  follows: 
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minimum 


where 


(Rj$/ft)  •  L 
Kf  •  (L/K) 1,0489 


0-4) 


RN/ft  is  determined  from  standard  atmospheric  tables  or 
is  input. 

L  is  the  characteristic  length  of  the  component. 

K  is  the  admissible  surface  roughness  and  is  an  input- 
quantity. 


and, 

«  37.587  +  4.615M  +  2.949M2  +  4.132M3 


For  mixed  laminar- turbulent  flow,  transition  location  is 
specified  for  the  upper  and  lower  surfaces  of  the  wing.  For  the 
laminar  portion  of  the  flow,  the  Blaoius  skin- friction  relation 

Cf  -  1.328//R^  ■  (3-5) 

1  Laminar 


n 

v  *• 


-4-0  19 


-•  W*  *.  J 


-.12 


■J  e 


VAU  VU  U  C»  CUW  LJ.QU0  1.' 


where  Cf/Cf^ 

tion  point.  At  the  transition  point,  Xr»  the  laminar  momentum 
thickness  is  matched  by  an  iterative  process  to  a  turbulent 
momentum  thickness,  which  begins  some  fictitious  distance,  AX, 
ahead  of  transition.  The  skin- friction  coefficient  for  the 
turbulent  part  of  the  flow  is  calculated  from  Equation  3-3, 
where  the  Reynolds  number  is  calculated  from 


RNl  *  (AX  +  L  -  Xr) • (RN/ft) 


TVio  *ia1i 


/*. 


r.rl  tV»  f-i-oin  *3  on 

W  J-  VA»  VA.UL1U  ^  C 


■P4  1 

■A.  1.11U  J.  Jf 


given  fcy 


Cf 


L-Xr. 

L  ^  CfTurb. 


(3-6) 


(3-7) 


Calculated  values  of  Cf  versus  Rjj  are  presented  in  Figures  2 
through  7  for  mixed  laminar- turbulent  flow. 


3.1.2  Form  Factors 


The  component  form  factors,  FF,  account  for  the  increased 
skin  friction  caused  by  the  supervelocities  of  the  flow  over 
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FIGURE  5  SKIN  FRICTION  ON  AN  ADIABATIC  FLAT  PLATE 


the  body  or  surface  and  the  boundary- layer  separation  at  the 
trailing  edge.  The  form  factor  for  the  "body"  component  is 
computed  as 

FF  «  1  +  60/FR3  +  0,0025  «  FR  (3-8) 


where 


Fr  =»  Component  Length 
y/ Width  x  Height 

For  "nacelle"  components,  the  form  factor  is  given  by 

FF  -  1  +  0.35/FR  (3-9) 

Equations  3-8  and  3-9  were  obtained  from  the  Convair  Aerospace 
Handbook  (Reference  4)  and  also  appear  in  the  DATCOM  (Reference  2). 

The  airfoil  form  factors  depend  upon  airfoil  type  and  stream- 
wise  thickness  ratio.  For  6-series  airfoils,  the  form  factor  is 
given  by 

FF  =  1  +  1.44(t/c)  +  2(t/c)2  (3-10) 

For  4-digit  airfoils,  the  form  factor  is  given  by 

FF  *  1  +  1.68(t/c)  +  3(t/c)2  (3-11) 

For  biconvex  airfoils,  the  form  factor  is  given  by 

FF  -  1  +  1.2(t/c)  +  100(t/c)4  (3-12) 

And  for  supercritical  airfoils,  the  form  factor  is  given  by 

FF  *  1  +  K]Cld  +  1.44(t/c)  +  2(t/c)2  (3-13) 

The  factor  KiC]^  in  Equation  3-13  is  an  empirical  relationship 
which  shifts  the  6-series  form-factor  equation  to  account  for 
the  increased  supervelocities  caused  by  the  supercritical- 
section  design  camber  Cj<j.  The  factor  Ki  (derived  f  ’om  experi¬ 
mental  data)  is  shown  plotted  in  Figure  8  as  a  function  of  the 
Mach  number  relative  to  the  wing  Mach  critical.  Equations  3-10 
and  3-11  were  obtained  from  informal  discussions  with  NASA/LRC 
personnel;  Equation  3-12  appears  in  both  the  DATCOM  and  the 
Convair  Handbook. 
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Wing  Compressl 


3.1.3  Interference  Factors 


Hie  component  interference  factors,  IF,  account  for  the 
mutual  interference  between  components.  For  the  fuselage,  the 
interference  factor  is  given  by 

IF  -  Rw_b  (3-14) 

where  Rw-B  *-8  shown  plotted  in  Figure  9  as  a  function  of  fuse¬ 
lage  Reynolds  number  and  Mach.  For  other  bodies  such  as  stores, 
canopies,  landing  gear  fairings,  and  engine  nacelles,  the  inter¬ 
ference  factor  would  be  an  input  factor  based  on  experimental 
experiences  with  similar  configurations.  The  Convair  Aerospace 
Handbook  (Reference  4)  recommends  using 

IF  *  1.0  for  nacelles  and  stores  mounted  out  of 
the  local  velocity  field  of  the  wing 

IF  ”  1.25  for  stores  mounted  symmetrically  on 
the  wing  tip 

IF  *  1.3  for  nacelles  and  stores  if  mounted  in 
moderate  proximity  of  the  wing 

IF  *  1.5  for  nacelles  and  stores  mounted  flush 
to  the  wing  or  fuselage. 

The  interference  factor  for  the  main  wing  is  computed  as 
IF  *  RLS  Rw-b  0-«> 

where  R^-b  the  wing-body  interference  factor  presented  in 
Equation  3-14,  and  Rls  1®  the  lifting  surface  interference 
factor  presented  in  Figure  .10.  For  supercritical  wings  the 
wing  interference  factor  is  set  equal  to  one.  Other  airfoil 
surfaces  such  as  horizontal  or  vertical  tails  use  an  interfer¬ 
ence  factor  determined  by 

IF  -  Rls  •  Hf  (3-16) 

where  Hf  is  the  hinge  factor  obtained  from  input  (use  Hf  »  1.0 
for  an  all  movable  surface,  1.1  if  the  surface  has  a  flap  for 
control).  The  factors  Ry-B  and  RLS  are  plotted  in  Reference  3 
and  also  appear  in  the  DATCOM. 


26 


3.2  Camber  Drag 

Hie  minimum  drag  contribution  of  the  wing  twist  and  camber 
is  related  to  the  lift  coefficient  of  the  polar  displacement, 
ACl>  by  the  equation 

C°CAMBER  "  T^e  (3-17) 


This  increment  is  called  camber  drag  and  represents  a  drag  incre¬ 
ment  between  minimum  profile  drag  and  Cd^n*  The  span  efficiency 
value,  e,  is  related  to  the  induced  drag  factor,  K,  by  the  equa¬ 
tion 

1 

6  "  7T-  AQ.  •  K 


If,  for  some  reason,  e  £  1,  an  alternate  equation,  obtained  from 
Reference  4,  is  used: 


CdCAMBER 


(3-18) 


3.3  Base  Drag 


Data  presented  in  Reference  5  were  used  to  establish  equa¬ 
tions  from  which  the  base  drag  of  bodies  could  be  determined. 

The  trends  of  these  data  show  three  different  phases:  (1)  a 
gradual  rise  of  CDjjase  at  transonic  speeds  up  to  M  *  1,  (2)  a 
relatively  constant  drag  level  supersonically  up  to  about  M  * 
1.8,  and  (3)  a  steadily  decreasing  value  of  drag  above  M  »  1.8. 
The  resulting  empirical  equations  are  given  as 


f  S 

(0.1  +  0.1222M8)  ~-ae,M  i  1 

sRef 

CDBase  “  <  0.2222SBase/SRef  ,1.0<M  <  1.8  (3-19) 

(1-42sBase/SRef)/(3-15  +  M2),  m  *  l-8 
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3 , 4  Wave  Draj 


Supersonic  wave  drag  is  determined  on  the  basis  of  a  compo¬ 
nent  buildup  for  which  simplified  shapes  are  assumed.  Three 
basic  simplified  shapes  are  used  to  represent  the  airplane: 
bodies,  nacelles,  and  wings.  The  component  buildup  assumes 
that  the  total  wave  drag  is  the  sum  of  the  isolated  wave  drag 
of  each  component  and  does  not  allow  for  the  mutual  interfer¬ 
ence  between  components.  However,  the  component- buildup  method 
does  give  wave-drag  results  comparable  to  average  configurations 
which  have  some  degree  of  favorable  and  unfavorable  interference. 


3.4.1 


The  technique  used  to  estimate  wing  wave  drag  evolved  from 
a  method  that  applies  transonic  similarity  theory  to  straight 
wings.  Data  correlations  at  Mach  1  were  performed  on  a  large 
number  of  unswept  wing  configurations  with  blunt  and  sharp 
leading-edge  airfoils.  For  the  Large  Aircraft  program,  these 
results  were  represented  by  an  analytical  function  common  to 
both  types.  The  equations  were  then  modified  for  M  >  1  to  pro¬ 
duce  a  peak  value  at  low  supersonic  speeds  and  then  to  decrease 
at  high  Mach  numbers  to  values  predicted  by  straight-wing  linear 
theory  for  equivalent  two-dimensional  configurations.  Finally, 
sweep  effects  were  included.  The  resulting  semi-empirical  equa¬ 
tions  are  presented  below: 

~  _ 2KtKwKcKb _ _ 

D  3  ’  r  1/aS*  aV  1 

/JKbVCV"*  - 3T«CT  +  — - 3 


|  1+(1+A)§/ 


1+ 


where 


_ 3.33  KtKfr-KcKb _ 

w«3-8<v“  + 

1+(|  +k)Fp0 


cD/(t/cr'J 

airfoil  thickness  distribution  factor 

l  +  4[i_|i(l+iV|2>]2-±V^  (1- 


1+3AR, 


(3-20) 


h)2 

c  ' 
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Kfo,  airfoil  factors 

"  1.0,  *  1.2  for  double-wedge  sections 

K5  ■  1.069,  K*,  ■  1.0  for  curved- type  sections 

Xt  ™  location  of  section  maximum  y  ordinate 

rQ  *  section  leading-edge  radius 

Kc  "  airfoil  camber  factor 

1  +  f(h/t)2 

h  *  section  camber  (maximum  y  ordinate) 


Kp 


1  /  *y 

<^osAle  +(YTA)?(tan2A^  "  tan2ATE> 
1  +  - 


(1+A2) 

0.3  +  0.7  Kp 


1  2 
2(tanALE  +  tanAfg) 


(1+2A^) 


9 


^/(t/c)1/3  -  >/i?Ti/(t/c)1/3 


^  >  ^lim 
&  ^  ^lim 


cosA^g  +-  cosArpE 


^lim  "  |taiu‘vL E 


and  where  ARe  is  the  straight-wing  AR^ having  the  s an 
Cj)/(t/c)3'3  at  M  -  1.0  as  Ak,  where  AR  ■  AR(t/c)^'  . 


me  value  of 


A*  , 

The  value  of  ARe  is  determined  by  solving  the  following 
equation  by  use  of  an  iterative  method: 

2  +  3.33  _  /  2  +  3.33 

4^  +  ARe3  4-5  +  1  V4-  +  AR3  4t+1 

ARe  AReJ  AR  AR-3 
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The  term  fi/im  represents  the  approximate  value  of  ft  ai:  which 
CD/(t/c)->/J  will  maximize,  provided  the  body  is  essentially 
cylindrical  where  the  wing  is  attached.  If  the  body  is  area- 
ruled,  the  peak  value  of  Cp/(t/c)5'3  may  or  may  not  be  closely 
approximated. 


3.4.2  Body  Wave  Drag 

The  fuselage  body  wave  drag  is  computed  by  dividing  the 
body  into  two  parts,  consisting  of  a  simplified  pointed  nose  and 
a  simplified  boattail.  That  is. 


Ci>wbody 


amax 

SREF 


+  °dPBT 


amax 

SREF 


(3-21) 


Nose, wave  drag,  Cpp  ,  is  determined  from  Linnell's  empirical 
equation  K 


(£„2+l)CD 


Pit 
*  n 


1.2  + 

1  +  1.9(^/ Vf^+l) 


(3-22) 


for  the  supersonic  wave  drag  of  parabolic  noses  (Reference  6). 

For  Mach  numbers  between  1.2  and  1.0,  the  nose  wave  drag  is 
determined  from  the  curves  of  Figure  11,  which  were  derived  from 
the  transonic  drag  rise  of  ogive  noses,  as  presented  in  Figure 
III. B. 10-9  of  Reference  4,  and  using  Equation  3-22  as  a  supersonic 
limit.  The  nose  fineness  ratio,  fN,  is  calculated  from  the  nose 
length,  and  the  maximum  cross-sectional  area,  A^^x,  as 

fN  ^^W^Amx 


Boattail  wave  drag,  Cp™,  is  determined  as  a  function  of 
boattail  fineness  ratio  (fg),  base  diameter  to  maximum  diameter 
(dg/d),  and  Mach  number.  This  is  done  by  computing  CppBT  at 
five  values  of  (dg/d)  and  interpolating  to  the  desired  value. 
The  general  form  of  these  equations  is  given  below: 
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For  /fN  ^  1  and  dB/d  -  •  ~  , 

Cdpbt(D  -  -^|  ao(I)+A  (I)-/9/fB+A2(I)-(/3/fB)2+A3(I)-(p/f3)3 

fB'  L 

(3-23) 

For  ^/fN  >  1  , 

CdPBT(I)  ™  ~2  A4(I)-(fN/^)  (3-24) 

fB 


The  polynomlnal  coefficients  of  Equations  3-23  and  3-4*4^ 
determined  from  a  least-square  fit  of  Fig.  III. B. 10-9  of  Refer¬ 
ence  4  for  ogive  boattails  ,  are  tabulated  below: 


I 

dB/ d 

Ao 

Al 

a2 

a3 

1 

0 

1.165 

-0.5112 

-0.5372 

0.3964 

0.513 

2 

0.4 

1.067 

-1.709 

1.6632 

-0.686 

0.3352 

3 

0.6 

0.7346 

-1.4618 

1.5795 

-0.6542 

0.198 

4 

n  a 

n  osss 

_n  snrm 

v/  »  ^  v/  v/  vy 

n  sno/, 

v  •  W  *.-V 

_n  0077 

v  #  fc  v  r  # 

ft  ft/.Q/. 

V  •  VT/-T 

5 

1.0 

0 

0 

0 

0 

1 

O 

3.4.3  Nacelle  Wave  Drag 

The  nacelle  wave  drag  is  calculated  by  a  method  similar  to 
that  used  for  the  fuselage: 


Nac 


(Cjw  +  CD_J  ^ 

\jn  foi  wKEi*’ 


(3-25) 


The  equation  used  to  calculate  for  open- nose  bodies  is 

%N  "  [U  -  rIH/R>/%]1‘5/^ 


(3-26) 


where 


fir 


r  INLET  "  ^AINLET 

R  “  'fAMAX^^* 

This  equation  is  a  curve  fit  of  Figure  III. B. 10-6  of  Reference  4. 
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3.5  Fuselage  Aft- End  Uplweep  Drag 

The  main  parameters  affecting  the  fuselage  aft-end  upsweep 
drag  are  the  upsweep  angle  and  the  crossflow  drag  coefficient 
of  the  rear  fuselage  sections  in  the  local  flow,  including  modi¬ 
fication  by  wing  downwash.  Data  in  Reference  7  indicate  that  the 
afterbody  drag  increases  rapidly  with  upsweep  angle  B  ,  but  decreases 
with  increasing  fuselage  angle  of  attack.  The  curves  in  Figure  12, 
obtained  from  the  data  presented  in  Reference  7,  are  used  to  predict 
aft-end  drag  as  a  function  of  angle  of  attack. 


3.6  Miscellaneous  Drag  Items 

In  the  preliminary  design  stage  of  aircraft  drag  estimation, 
the  drag  due  to  surface  irregularities  such  as  gaps  and  mismatches, 
fasteners,  small  protuberances,  and  leakage  due  to  pressurization 
are  estimated  by  adding  a  miscellaneous  drag  increment  which  would 
be  some  percentage  of  the  total  friction,  form, and  interference 
drag*.  The  miscellaneous  drag  varies  between  10  and  20  percent  of 
the  total  friction,  form,  and  interference  drags  for  typical  air¬ 
craft.  The  Large  Aircraft  program  computes  miscellaneous  drag  by 
use  of  a  percentage  factor  specified  as  input  to  the  program. 
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4.  DRAG  DUE  TO  LIFT 


The  Large  Aircraft  program  predicts  drag  due  to  lift  by  one 
of  several  methods,  depending  on  the  aerodynamic  conditions  at 
which  a  solution  is  desired.  The  various  regions  are  illustrated 
in  Figure  13;  they  are  discussed  in  the  following  subsections  in 
the  numerical  order  shown  in  the  figure. 


McRq  ML1  **L2 
MACH  NUMBER 


Figure  13  Lift  and  Speed  Regions  for  Calculation  of 
Drag  Due  to  Lift 


4.1  Subsonic  Polar  Prediction  below  Polar  Break 


Region  1  is  bounded  by  the  critical  Mach  number  and  by  the 
Cl  at  which  the  polar  break  occurs  (CLb»).  In  this  subsonic, 
low-lift  region,  the  drag  due  to  lift  can  be  determined  from 


(4-1) 


Cdl  “  K(CL  -  acl)2 

where  the  drag  due  to  lift  factor,  K,  is  predicted  from 


1-R  R 

Cl«  ^AR  e0 


(4-2) 


In  this  equation,  a  leading-edge  suction  parameter  R,  is  used  to 
relate  K  to  the  lower  bound  of  drag,  1/ 3TAR,  for  full  leading- 
edge  suction  (R-1.0)  and  to  the  upper  bound  of  drag,  1/Cl„,  for 
zero  leading-edge  suction.  Body  effects  are  accounted  for  in 
Equation  4-2  by  computing  eQ,  shown  plotted  in  Figure  14  as  a 
function  of  taper  ratio  and  body-diameter- to-span  ratio  (d/b). 


The  correlation  of  leading-edge  suction  on  induced  drag  was 
first  developed  by  Frost  (Reference  8)  and  was  later  extended 
for  additional  planform  effects  and  higher  subsonic  Mach  numbers 
(Reference  9).  A  study  by  NASA  (Reference  10)  showed  that  air¬ 
foil  camber,  conical  camber,  sharp  leading  edges,  leading-edge 
flaps,  Reynolds  number,  and  sweep  have  significant  effects  on 
the  suction  parameter.  H.  John  (Reference  11)  improved  the 
correlation  of  R  for  plane  wings  at  low  Reynolds  number  by  in¬ 
cluding  airfoil  thickness  along  with  leading-edge  radius. 


The  procedure  followed  in  the  Large  Aircraft  program  to 
determine  R  is  as  follows: 


1.  Using  the  leading-edge  radius  and  the  leading-edge 
sweep  for  each  wing  panel,  determine  f)  from  the 
equation 

rNLEr  x10"'J  co <Ale  V 1- M2cos ;  ALE  >  20° 

_  (4-3) 

rNLEr  x10"3(5-6.51lALE)/l^Mrcos^.LE;  Ale  <  20° 

(Ale  *-n  radians) 

The  switch  from  the  cotangent  term  is  made  to  prevent 
flfrom  going  to  infinity  as  sweep  approaches  zero. 

The  value  of  fi  is  then  used  to  read  Rj  from  Figure 
15,  which  is  a  plot  of  leading-edge  suction  for  thin, 
round-nose,  uncarabered  airfoils  developed  in  Refer¬ 
ence  9. 
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Effect  of  Reynolds  Number  on  Leading-Edge 
Suction  for  Blunt -Uncambered  Airfoils 


2.  Using  the  leading-edge  sweep  of  the  wing  panel, 
determine  RmIH  fro®  Figure  16.  The  plot  of 
Figure  16  was  obtained  from  the  results  reported 
in  Reference  10  for  sharp -leading -edge  wings.  The 
leading-edge -suction  value  for  sharp- leading-edge 
wings  is  Independent  of  Reynolds  and  Hact|  numbers . 

3.  If  the  value  of  Rj  determined  from  step  (1)  Is  less 
than  the  value  of  Rmir  determined  from  step  (2),  set 
Rj  equal  to  Rmin- 

4.  Using  the  wing  panel  thickness  and  leading-edge 
radius  Reynolds  number,  determine  a  thickness  cor¬ 
rection  to  leading-edge  suction,  R-j,  from  Figure 

17.  Figure  17  was  determined  from  the  data  presented 
in  Reference  11.  The  increment  in  suction  parameter 
i6  then  added  to  the  value  of  Rf  determined  from 
step  (3). 

5.  Determine  the  effect  of  either  section  camber  or 
conical  wing  camber  on  the  R  factor  from 


R-i  -  R-r  +  f0.824-RT)  •  (Ci  .+€?  ) / 0 . 6 

*  j  '  a'  x  uCOn 


tt.  /.  \ 


If  Ri  >  0.874,  Ri  -  0.874. 

Correlation  with  experimental  data  and  the  results 
of  Reference  10  indicate  that  R  does  not  decrease 
as  much  for  low  Reynolds  number  when  the  wing  is 
cambered  compared  to  an  uncambered  wing. 

The  accompanying  sketch  ^ 

shows  the  relative  effect  ^ 

that  Equation  4-4  will 
predict  for  cambered 
wings . 


**LER 

6.  Obtain  the  effective  R  for  the  composite  wing 
from  a  span-weighted  average  of  the  individual 
for  each  panel  as  follows : 


Ref.  10,  Fig.  4 


Figure  16  Leading-Edge  Suction  for  Sharp  Airfoils 
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Figure  1?  Effect  of  Reynolds  Number  and  Thickness  on 
the  Leading-Edge  Suction  Factor 
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i-  z 


Ri 


YW(l-fl)  -  Y W(t) 

(b/2  -  d/2) 

% 


where  W(i+1)  -  W(i)  is  the  width  of  the  ith  panel. 

7.  Using  Figure  18  and  the  parameter  AR  *  ^ /cosAle*  calcu¬ 
late  the  effect  of  wing  planforro  on  R.  The  increment  AR 
due  to  planform  is  then  added  to  IT,  determined  from  step 
(6),  to  obtain  the  final  R  value  used  in  Equation  4-2. 

The  polar  displacement,  A  Cl,  Is  related  to  the  lift  coeffi¬ 
cient  for  minimum  profile  drag,  CLgpj»  by  the  equation 


aCl  -  (i  -  ttaTk)Clopt 

The  lift  coefficient  for  minimum  profile  drag  is  related  to  the 
camber,  twist,  and  asymmetry  of  the  configuration.  Figures  19 
and  20  present  data  (Reference  13)  for  the  effect  of  NACA  camber 
and  conical  camber  on  Clopt.  For  supercritical  wings,  the 
limited  amount  of  data  available  for  correlation  indicate  that 


0.75 

%PT  =  °*5195<CLd) 

where  Cl^  is  the  wing  section  design  lift  coefficient. 


(4-5) 


4.2  Supersonic  Polar  Prediction  below  Polar  Break 


The  drag  polar  in  the  supersonic  region  beyond  the  second 
limit  Mach  number  below  polar  break  (Region  2)  is  predicted  by 
Equation  4-1,  where 


1-R  +  % 

Cue  ?TAR  e0 


(4-6) 


This  equation  is  similar  to  Equation  4-2  for  the  subsonic  induced- 
drag  factor  except  for  the  use  of  R,  which  is  a  transonic  leading- 
edge  suction  factor.  For  Mach  numbers  greater  than  Mach  critical, 
the  suction  factor  is  predicted  from 


R  -  R0/(l  +  nAM  +  (nAM)2) 


(4-7) 
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Ref.  13 


U-c- 


FLgure  19  Lift  Coefficient  for  Minimum  Profile  Drag 
NACA  Camber 

DESIGN  MACH  -  1,0 
0.8  b/2  RAY  LINE 
•6  Ref.  13 


.3  .4  -5 

Conical  Camber  -  Cl^q^ 


Figure  20  Lift  Coefficient  for  Minimum  Profile  Drag 
Conical  Camber 
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where  Rq  ■  leading-edge  suction  factor  at  the  critical 
Mach  number  at  zero  lift, 

n  »  12(coaALE)1,6 

AM  -  M  - 

The  variation  of  (R/Rq)  for  Mach  numbers  greater  than  Mq^  is 
shown  in  Figure  21.  This  method  of  predicting  polar  shape  factor 
produces  a  continuous  decrease  in  leading-edge  suction  so  that 
in  the  limit;  as  the  Mach  number  approaches  the  sonic  leading-edge 
condition,  the  polar  shape  approaches  (1/Cl^. 

The  supersonic  polar  displacement  for  NACA  camber  is  calcu¬ 
lated  from 


^CLd(0.25-0.225^otjtLE);  fcot^LE  4:  1.11 
AcL  “  ]  .  (4-8) 

(  0;  fcotj^g  V  1.11 

and  the  supersonic  polar  displacement  for  conical  camber  is 
calculated  from 

(cLc^o.m-o.i^^);  fcocAie  l.n 

ACf]  (4-9) 

j0;  pcotJjj.^1.11 


These  equations  were  obtained  from  a  simple  curve  fit  of  the 
data  presented  in  Figures  94  through  97  of  Reference  13. 


4.3  Transonic  Polar  Prediction 

In  the  transonic  region  bounded  by  Mach  critical  (M^g)  and 
the  first  limit  Mach  number  (MEl)  (Region  3),  the  induced  drag 
is  computed  by  adding  drag  rise  to  the  basic  polar: 

cDi.  "  K<cl  -  4cL>2  +  CDrcl  (4-10) 

The  basic  polar  shape  is  calculated  up  to  Mq^  in  the  same  manner 
as  described  for  Region  1.  Beyond  the  basic  polar  does  not 

change.  An  incremental  drag-rise  term  )  is  calculated  as  a 
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function  of  lift  and  is  added  to  the  basic  polar  to  determine 
the  total  drag  due  to  lift.  A  complete  description  of  the 
techniques  used  to  determine  the  drag-rise  increment  is  given 
in  Subsection  5.3. 

In  the  transonic  region  above  the  first  limit  Mach  number. 
Mm  (Region  4),  the  drag-rise  term  in  Equation  4-10  becomes  less 
accurate.  Therefore,  the  drag  polar  in  Region  4  is  calculated 
by  interpolation  between  the  polar  shape  factors  calculated  in 
Regions  3  and  2.  The  equations  for  K  and^Ci,  are  given  by 


M-Ml1 

K  -  KU  +  (KL2  -  KU)  Ml2_„li  (4-11) 

M-M,, 

ACL  ■■  ACLu  +  (ACLl2  -  ACLu)  (4-12) 

where  Mji  <  M  <  Mro-  "T*16  P°lar  shape  factors,  and  Clu» 
are  determined  from  a  least-square  curve  fit  of  the  polar  shape 
computed  by  Equation  4-10  at  Mu*  The  polar  shape  factors  Kl2 
and  are  computed  from  Equations  4-6  and  4-8  at  Mu- 

The  limit  Mach  numbers  are  determined  from 

ML1  *  +0-05 

If  MU  <  0.95,  then  Mu  -  0.95 

If  ML1  >  1.00,  then  ML1  -1,0 

mL2  =  MLi  +0.15 


4 . 4  Subsonic  Polar  Prediction  above  Polar  Break 

The  polar  region  between  the  polar-break  lift  coefficient 
(CLpg)  and  the  initial-stall  lift  coefficient  (C^g)  is  the 
region  in  which  leading- edge  separation  and  reattac'nment  occurs, 
causing  the  polar  to  deviate  from  a  parabolic  shape  (Region  5 
in  Figure  13).  Whether  or  not  this  region  exists  (i.e.,  the  flow 
reattaches  after  separation  and  allows  the  wing  to  reach  a  higher 
lift  coefficient  before  final  separation  occurs  at  the  trailing 
edge)  is  de  ;ermine 3  by  the  type  of  airfoil,  the  Reynolds  number,  and 
the  leading-edge  wing  sweep.  For  thin  wings,  low  Reynolds  numbers, 
or  highly  swept  wings,  the  values  of  CjjT>g  and  CLpB  are  et?ual* 


Leading-edge  sharpness  is  a  measure  of  the  type  of  separation 
likely  to  occur.  Blunt,  thick  airfoils  generally  exhibit  trailing- 
edge  separation,  while  very  thin  airfoils  exhibit  leading-edge 
separation.  Airfoils  of  moderate  thickness  are  likely  to  separate 
and  reattach  at  the  leading  edge,  followed  by  trailing-edge  sepa¬ 
ration  (stall)  at  higher  lift  coefficients.  Associated  with  the 
leading-edge  separation  and  reattachment  is  a  loss  in  leading-edge 
suction,  which  produces  an  increase  in  drag  due  to  lift.  Above 
Cldb»  t*ie  separates  completely  along  the  wing  and  the  drag 

increases  more  rapidly. 

The  prediction  method  in  the  Large  Aircraft  program  utilizes 
the  sharpness  parameter  of  the  airfoil,  Ay,  as  defined  in  Equation 
2-28.  If  Ay  is  less  than  or  equal  to  1.65,  leading-edge  flow 
separation  is  assumed  to  occur.  Also,  if  the  leading-edge  sweep 
is  greater  than  or  equal  to  50  degrees,  it  is  assumed  that  lead¬ 
ing-edge  separation  occurs.  If  Ay  is  greater  than  2.05,  a  lead¬ 
ing-edge  separation  and  reattachment  occurs,  followed  by  a  trail¬ 
ing-edge  separation.  For  values  of  Ay  between  1.65  and  2.05,  a 
transition  region  exists  in  which  the  behavior  varies  between 
the  condition  of  full  leading-edge  separation  at  Ay  “  1.65  and 
full  leading-edge  flow  reattachment  at  Ay  «  2.05. 

Because  the  polar-break  lift  coefficient  is  a  function  of 
many  variable?,  it  has  proved  to  be  a  difficult  quantity  to 
predict.  Data  correlations  at  subsonic  and  transonic  speeds 
performed  during  the  development  of  the  Large  Aircraft  program 
indicate  that  the  polar-break  point  can  better  be  determined 
with  angle  of  attack  as  the  parameter  rather  than  lift  coeffi¬ 
cient.  These  correlations  resulted  in  a  method  that  determines 
the  angle  of  attack  at  polar  break  as  a  function  of  Mach  number. 

Ay,  sweep  angle,  and  wing  camber.  Consequently,  the  polar  break. 
Cl,  is  calculated  as 

CLpb  =  CLa  *  (aPB0  +AQPb)  (4~  13) 

where  (Gpgo/co*^c/4^  Is  shown  plotted  in  Figure  22  as  a  function 
of  Ay  and  M  cosAc/4*  The  term  AGPB  accounts  for  section  camber 
ani  is  determined  from 

AOpb  -  (12.05-4.1  M  ““4/4) (cosj°  ) 

which  was  derived  principally  from  correlating  the  experimental 
data  in  Reference  12.  lor  wings  with  conical  camber,  an  Incre¬ 
ment  in  ClP£  is  obtained  from  Figure  23. 
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Airfoil  Lsadtng-Edge 


(4-14) 


The  drag  due  to  lift  Is  expressed  In  Region  5  as 
CDl  -  K(CL-ACL)2  +  K'(CL-CLpB)2 

cLpfi  <  CL  ^  cLdb 

where  K'  «  0.518/  >/  AR^  (Reference  1). 

The  upper  boundary  of  Region  5  represents  the  lift  coeffi¬ 
cient  at  which  trailing-edge  separation  occurs.  It  is  predicted 
aa 

Ct-DB  -  CWb  +  1  Kc  •  O  -  aCL] 

where 

0  ;  Ay  ^  1*65 

T  *  (Ay  -  1.65)/0.4;  1.65  <AY<  2.05 

1.0  ;  Ay  ^  2.05 

and  (0Cl/5Rn)  and  0Cl  are  shown  plotted  in  Figures  24  and  25. 

The  drag  polar  above  Cldb  (Region  6  in  Figure  13)  increases 
sharply  from  the  subsonic  attached-flow  condition.  The  polar 
prediction  for  lift  coefficients  above  Cldb  *s  determined  by  a 
modification  of  the  empirical  method  presented  in  Reference  13 
whereby 

CDl  *  cDdb  +  kDcL2  +  acDB  (4-15) 

CL  >  cL])b 

where  Cddb  is  the  predicted  lift  dependent  portion  of  profile 
drag  at  Cldb»  is  the  theoretical  induced  drag,  andACpe 

is  a  correlated  separation  drag  increment  obtained  as  a  function 
of  CldB'  The  lift  dependent  profile  drag  at  Cldb  is  given  by 

Cddb  -  K(CLDB-AOL)2+K'(CLDB-CLpB)2-KDCLDB2  (4-16) 


where  the  theoretical  induced  drag  factor,  Kp,  is  predicted  from 


n  AR  e 


(4-17) 


where  e'  is  a  modification  of  the  classical  theoretical 
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drag-due- to- lift  factor  1/ 7TAR  to  account  for  nonelliptical  span 
loadings  and  body  effects.  The  factor  e~  is  calculated  from 

e'  ■  e'w  t1  -  <d/b>2]  (4-18) 

where  the  wing  planform  efficiency  factor  e'w  is  as  plotted  in 
Figure  26  as  a  function  of  taper  ratio,  sweep,  and  aspect  ratio. 
The  data  of  Figure  26  were  obtained  from  a  Weissinger  lifting- 
line  solution  presented  in  Reference  13. 

The  drag  above  Cl™  represents  the  separation  drag  component 
when  major  separation  effects  are  present.  Simon  et  al.  (Refer¬ 
ence  13)  measured  this  drag  relative  to  the  profile  drag  at 
the  drag-break  lift  coefficient  and  present  correlated  curves 
of  ACdb  versus  Cl>  Cl™,  end  Mach  number.  The  ACddb  data  were 
curve  fitted  and  result  in  the  equation 

ACdb  -  Kb(CL-CLdb)2+.08  AB(CL-CLi)B)2  ;  CL  >  CLob  (4-19) 

The  factor  Kg  is  shown  plotted  in  Figure  27.  The  program  does 
not  vary  Kg  with  Mach  number  since  the  drag-rise  term  is  included 
in  the  polar  buildup  above  Cldb* 

For  conventional  wings  the  polar  predicted  by  Equation  4-15 
is  continued  up  through  Clmax*  F°r  ^.ow-aspect-ratio  or  cranked 
wings  (see  Section  6)  that  develop  vortex  lift  at  the  higher 
lift  coefficients,  the  zero-suction  drag  polar  predicted  by 

cDl  "  CL  tana  (4-20) 

is  compared  against  the  drag  polar  predicted  by  Equation  4-15 
and  the  drag  due  to  lift  is  set  equal  to  whichever  is  lowest. 

This  produces  a  drag  polar  as  shown  in  the  sketch  below.  The 
Cl- a  variation  for  vortex  lift  is  predicted  by  the  methods 
given  in  Section  6. 
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Cl*  tana 


4 . 5  Supersonic  Polar  Prediction  Above  Polar  Break 

Supersonic  Region  7  above  polar  break  is  predicted  by  the 
semi-empirical  method  developed  in  Reference  13.  The  equation 
for  predicting  supersonic  drag  due  to  lift  above  polar  break  is 
given  by 

Cdl  =  (K-K* )(Clpb-  ACl)2+K' (Cl-  *Cl>2;  CL  >  CLpB  (4-21) 

where  K  and  A Cl  ate  the  polar  parameters  in  the  low- lift  region 
discussed  in  Subsection  4.2. 

The  polar-break  lift  coefficient,  CLpg*  is  correlated  as  a 
function  of  sweep,  aspect  ratio,  camber,  etc.  in  Reference  13. 

A  curve  fit  of  the  data  in  Reference  13  results  in  the  following 
equations  for  CLpg- 

cLpg  *  CLS1+1.25(CLS9+^CLS9-CLSi)(?cotJ^LE"*1)'H:,,5CI'd  (4-22) 

where  the  factors  CLSi»  Cls9>  and  ACls9  are  shown  plotted  in  Figure 
28  as  a  function  of  aspect  ratio. 
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The  polar  shape  factor,  K' ,  above  C^pg  is  computed  in 
Reference  13  as 


K' 


where 


H 


H 


1.1;  ARtan\LE  ^3.5 

1.1  +  0.1(ARtanALE  “  3.5);  AR  taii^E  ^3.5 


(4-23) 


(4-24) 


59 


5.  CRITICAL  MACH  NUMBER 
AND  DRAG  RISE 


The  drag-divergence  Mach  number  or  Mach  critical,  is  defined 
as  the  Mach  number  at  which  a  rapid  drag  rise  intercepts  the  sub¬ 
sonic  trend  in  drag.  The  British  method  of  predicting  the  critical 
Mach  number  for  two-dimensional  airfoils  (Reference  1^  appears  to 
be  the  most  accurate  empirical  method  available.  The  British 
method  uses  the  Sinnott  "crest  criteria" .where  the  low-speed 
pressure  at  the  airfoil  crest  is  related  to  the  drag-divergence 
Mach  number. 

The  Large  Aircraft  program  uses  a  method  analogous  to  the 
British  two-dimensional  critical  Mach  number  prediction  procedure 
in  order  to  predict  the  critical  Mich  number  for  a  finite-aspect- 
ratio  swept  wing.  The  critical  Mach  number  is  defined  as  the 
value  o?  freestream  Mach  number  which  produces  a  local  supersonic 
flow  measured  normal  to  the  sweep  of  the  isobar  at  the  crest. 

The  local  Mach  number  normal  to  the  crest  isobar  is  taken  to  be 
1.02  for  conventional  airfoils  and  1.05  for  supercritical  airfoils 
in  order  to  define  freestream  critical  Mach  number.  The  sonic 
condition  at  the  crest  can  be  predicted  by  means  of  a  simple 
equation  in  which  the  incompressible  pressure  at  the  crest  of 
the  airfoil  and  compressibility  factors  are  used. 

The  value  of  1.02  local  Mach  number  for  the  weak  shock  at 
crest  condition  for  drag  rise  used  in  Sinnot’ ' s  transonic  air¬ 
foil  theory  (Reference  15)  was  established  empirically.  Refer¬ 
ence  14  shows  that  this  method  should  predict  Mq*  to  within 
+0.015  for  the  majority  of  conventional  two-dimensional  airfoils. 
However,  as  shown  in  Reference  14,  with  "peaky"  airfoils  (as  in 
the  supercritical  airfoil)  the  onset  of  rapid  drag  rise  may  be 
delayed  until  the  shock  is  substantially  downstream  of  the  crest. 
The  predicted  value  of  Mcr  based  on  a  local  Mach  of  1,02  at  the 
crest  may  thus  be  conservative  by  more  than  0.02,  and  a  local 
Mach  of  1.05  is  necessary  to  achieve  good  correlation. 

The  following  subsections  discuss  the  methods  used  to 
predict  the  pressure  distribution  around  an  airfoil  and  to 
determine  from  the  pressure  at  the  crest  along  with  the 
method  used  to  estimate  the  drag  rise  above  M^p. 
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5 • 1  Pressure  Coefficient  Calculations 

The  incompressible,  inviscid  pressure  distribution  around 
the  airfoil  is  first  defined,  from  which  the  pressure  at  the 
crest  can  be  determined.  The  method  of  Weber  (Reference  16) 
was  used  for  the  pressure  distribution  calculations.  This 
method  requires  the  airfoil  surface  coordinates  to  be  deter¬ 
mined  at  the  chordwise  locations  defined  by 


X(V)  ■  ^(1  +  cos  ^L)  where  0^  ^  N 


(5-1) 


N  may  be  any  integer,  but,  in  this  program,  N  is  set  equal  to  32. 

The  Weber  formula  is  essentially  a  second-order  linear 
theory  whereby  the  pressures  are  determined  from  multiplication 
of  the  matrix  of  thickness  and  camber  ordinates  of  the  airfoil 
by  a  matrix  of  constants  given  in  Reference  16. 

The  formula  for  the  incompressible  pressure  distribution 
on  an  infinite  sheared  wing  was  obtained  from  the  incompressible 
form  of  Equation  93  in  Reference  16,  resulting  in 


1  -  Ct 


S(2)(X)  +  S(5)(X)  1 

cos  _A  J 


xjjcosoC  1  +  S'  ^(X)cos_A  +  (X)cos_yU 


+  sin  *  cosAji  +  ? 

—  <_  CUS./L  a  /  j 

+jcoso«L  (X)siny^  +  (X)sinTlJ 

+  sin  oi  sinA{l  +  ^ 

+■  sin2_Acos2<>/ ]  1 - r  - - rr? - = 

1  rc(2)/v,+eV5)/vN1 


+  sin  04  sin. 


;(X)±SV  ;(X 
cosy^ 


For  Cp  or.  the  upper  surface  the  +  is  used,  and  for  the  lower- 
surface  Cp  the  -  is  used.  Also> 


s(1)W  ~ZS% 

Zt/L 

s(2)(x)  = 

/1=1 

Zt/W- 

S(3)(X) 

//%  N-l  \ 

s(4'(x)  *  2s  ^ 

>i*1 

^s/t 

s(5bxi  -  Vs<5L,> 

n  /  y  <  / — 

-S/-C 

>i~l 

Tables  of  the  Sv  2,^  matrices  of  constants  are  given  in  Reference 
16  ,  and  ZJ/W_  are  the  thickness  and  camber  distribution 

at  the  control  p^rnt/t  given  b 7 


1 

zt >t  =  +  y jl  ) 


and 


Z 


y  jl  ) 


where  and  are  the  upper-  and  lower-surface  ordinates 
defined  at  the  control  point  /U.  given  by 

Y(/l)  “  j  O  +  cos^> 
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5 . 2  Critical  Mach  Number  Calculation  from 
Crestline  Pressure 


The  procedure  followed  to  determine  for  swept  wings  is 
similar  to  the  procedure  outlined  in  Reference  14  to  predict 
the  for  two-dimensional  sections.  In  the  procedure,  Equa¬ 
tion  5-2  is  used  to  compute  the  pressure  distribution  around 
the  airfoil  at  a  sweep  angle  determined  from 


A  = 


arcos(cosAc/2)ni  Ac/2  ~  40° 

f  (cos/\c/2)n+.76604n 


arcos 


3  L 


1- 


A 


c/2 


>  40 


(5-3) 


whereAc/2  Is  the  wing  mid-chord  sweep  at  the  semi-span  of  the 
wing  and  the  factor  n  is  determined  from 


n 


AR 

1 .4+AR 


The  sweep  angle  A  represents  an  effective  isobar  sweep  at  the 
mid-span  region  of  the  wing  as  affected  by  the  root  and  tip 
regions  of  the  wing.  The  procedure  used  to  determine  M^r 
based  on  the  crest  pressures  is  as  follows: 


1.  Determine  a  chordwise  incompressible,  inviscid 
pressure  distribution  for  an  angle  of  incidence 
(a).  Integrate  the  pressure  distribution  to 
obtain  the  lift  coefficient  (Cl^) . 

2.  Determine  the  chordwise  position  of  the  crest 
for  each  a,  the  crest  being  defined  as  the 
point  at  which  the  airfoil  surface  is  tangential 
to  the  undisturbed  freestream  direction  (0  =0*). 

3.  Determine  the  incompressible  pressure  coeffi¬ 
cient  at  the  crest  (Cpcrest)' 

4.  Use  Cp  ^  to  determine  Mpp  from  the  relation 

rcrest 

(P/Pl)(1+0.2Mcr2coc^)3,5  _ l 
CPcrest  - " 

o .  7McR2/  %  i-McR2cos2A 

where  (P/Pt)  is  the  ratio  of  local  static  pres¬ 
sure  to  freestream  stagnation  pressure  as 
determined  from 
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r  -i3-5 

P/Pt  -  [H  0.2M/J 

where  H|  is  the  local  Mach  number  normal  to  the 
isobar  sweep  A  at  the  crest  of  the  airfoil.  Mj 
is  set  equal  to  1.02  for  conventional  airfoils 
and  1.05  for  "peaky"  or  supercritical  airfoils. 
Equation  5-4  uses  a  Prandtl-Glauerc  compressi¬ 
bility  factor  to  correct  the  incompressible 
pressure  coefficient  for  Mach  number  rather 
than  the  Karman-Tsien  factor  used  in  Reference 
14.  References  17  and  18  recommend  using  the 
Prandtl-Glauert  factor  instead  of  the  Karman- 
Tsien  factor  in  the  Mcr  prediction  method  for 
highly  cambered  airfoils  or  general  airfoils  at 
high- lift  coefficients.  The  relationship  deter¬ 
mined  by  Equation  5-4  is  plotted  in  Figure  29. 

5.  Use  the  Prandtl-Glauert  compressibility  factor 
($p)  evaluated  at  M^r  to  obtain  the  lift  coef¬ 
ficient  Cld  from 

4’LD  “  C Li//*D 

6.  Repeat  Steps  1  through  5  for  a  set  of  incidences 
in  order  to  obtain  a  drag-rise  boundary  from  the 
set  of  points  (€Lj)J  Mcr)- 

The  critical  Mach  number  predicted  by  the  above  six  steps 
is  prevented  from  exceeding  the  critical  Mach  number  of  the 
fuselage  alone  (shown  plotted  in  Figure  30).  For  aircraft 
that  are  not  area-ruled,  where  the  isobars  are  allowed  to 
unsweep  at  the  wing-fuselage  juncture,  the  method  would  tend 
to  overpredict  Mcr  when  the  value  approaches  the  fuselage  Mcr. 
The  prediction-versus-test  Mcr  correlation  shown  in  Figure  31 
is  thus  applied  for  conventional-wing  predictions. 


5 . 3  Drag  Rise 

For  Mach  numbers  less  than  Mcr  the  drag  increases  slowly 
with  increasing  Mach  number.  This  drag  component  is  known  as 
compressible  drag,  or  drag  creep.  Methodology  for  estimating 
this  component  of  drag  tor  conventional  or  supercritical  wings 
was  included  in  the  subsonic  drag  buildup  in  Subsection  3.1. 
For  Mach  numbers  greater  than  Mqr>  drag  rise  begins  and  in¬ 
creases  rapidly  with  Mach.  Figure  32  illustrates  the  drag 
bookkeeping  system  followed  in  the  Large  Aircraft  program 
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SOURCE 

NACA  RM  L52B03 


Correlation  of  Critical  Mach  Number  for  Con 
vent ion a i  Wings 


whereby  beyond  Mach  1.0  the  drag  rise  and  the  interference  plus 
form  drag  are  replaced  by  wave  drag.  Hie  drag  rise  is  separated 
into  two  components,  drag  rise  due  to  lifting  surfaces  and  drag 
rise  due  to  all  other  components  on  the  aircraft.  The  drag  rise 
due  to  lifting  surfaces  is  represented  by 

CDRL  "  PL  (M-MCR)2  (5-5) 


where 


PL  «  25  '  (t/c  +  2f/c)*  (cos^/2 )3 

The  factor  is  a  function  of  the  wing  section  thickness,  t/c, 
maximum  ordinate  of  the  camber,  f/c,  and  the  midchord  sweep, 
■Ac/2-  The  total  drag  rise,  at  zero  lift,  due  to  lifting  and 
non- lifting  components  is  determined  by 

“  ^(^McRo)2  +  a3(M"MCR0^3  (5-6) 


where  a2  and  a3  are  defined  to  produce  a  continuous  zero -lift 
drag  curve  between  Mcr0  and  Mach  1.0.  The  drag  rise  is  curve- 
fitted  to  begin  at  Meg  with  zero  slope  and  end  at  Macb  -  1.0, 
matching  the  value  and  slope  of  the  wave -drag  curve.  The 
coefficients  a£  and  33  are  calculated  from 


3(CD.„  -  CD  >  -  <!  -  “CR  >Cd' 


a-j 


W1 


F&I 


Wl 


(1  '  ^CR0^C^1  ~  ^(CDwl  ~  CpF&l^ 
(L  " 


where  CDp&-,and  Cd^i  represent  the  Mach  1.0  value  of  the 

wave  drag,  form  plus  interference  drag,  and  the  slope  of  the 
wave  drag,  respectively.  The  non-lifting  contribution  to  the 
drag  rise  is  assumed  to  begin  at  MCRo  anc*  not  vary  with  lift; 
using  Equations  5-5  and  5-6,  an  equation  for  total  drag  rise 
at  any  lift  condition  can  be  determined  from 

%  "  PL(M-MCR)2  +  (a2-PL>(M-MCR0>2  +  a3(M'MCR0>3  (5“7> 
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The  change  in  M^r  with  lift  causes  the  subsonic  drag  polar 
to  increase  after  M^r  (8ee  Figure  32).  For  bookkeeping,  the 
drag  rise  is  separated  into  a  minimum  drag  contribution  and  a 
contribution  due  to  lift.  The  increment  to  the  minimum  drag  is 
determined  from  Equation  5-6,  and  the  increment  at  lift  is 
determined  by 

C%CL  -  %  -  CDRo  <5-8> 
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6.  LIFT 


The  untrimmed  lift  of  an  aircraft  can  be  represented  by 
the  equation 

CL  *  CL,*  (°*"  ^Lo) 

For  moderate  to  high  aspect  ratios  and  moderate  sweeps,  the 
lift  equation  is  linear  with  0c)  so  that  the  lift-curve  slope 
(Clk)  is  constant.  The  total  lift-curve  slope  of  the  aircraft 
is  given  by 

CLoL=  <CL.<)w_b  +  (CLJt  +  <clJb  (6-l) 

which  is  the  sum  of  the  wing  (including  body  carry-over  effects), 
horizontal  tail,  and  the  forward  portion  of  the  fuselage. 

The  following  subsections  describe  the  methods  by  which 
Cl^  andoCix)  are  calculated  along  with  the  method  of  calculating 
the  lift  in  the  nonlinear  range  up  to  Clj*^* 

6 • 1  Wing  Lift-Curve  Slope 

The  value  of  (Cl*)w_jj  is  predicted  by  the  use  of  several 
rather  involved  semi-empirical  equations.  These  equations  were 
developed  to  predict  wing  lift-curve  slope  as  a  continuous  ex¬ 
pression  in  the  subsonic,  transonic,  and  supersonic  regions. 

The  value  of  (Clv)w_b  is  expressed  as 


<CL, 


•^W-B 


<CL*> 


Basic 


Qr*T 

^rij/uN 

SrEF 


(6-2) 


where  (Ci^)Bas^c  is  the  wing-alone  Cl*  with  no  thickness  effects. 
The  factors  Kt  and  K^,  account  for  the  effect  of  airfoil- thickness 
plus  camber  and  fuselage  interference,  respectively.  The  equation 
for  (Cl*) Basic  was  evolved  from  the  Polhamus  (Reference  1  9)  equa¬ 
tion  for  trapezoidal  wing  . 


cl*  - 


(~)AR/5  7.3 


(6-3  ) 


2tiJ 


/  aO  2 

r  2] 

f(sr>+ 

l-(Mcoa/lc/r2)  i 

AR 
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for  subsonic  flow  and  the  linear-  theory  level  of 


cl* 


4/57.3 

'HTi 


(6-4) 


for  supersonic  flow  at  the  leading  edge  (M  >  1/ cosAy?)  • 

To  extend  the  Polhamus  equation  for  use  with  non-trapezo:  dal 
wing  planforms,  Spencer  (Reference  20)  replaces  cos^/2  with  the 
effective  cosine  mid-chord  sweep  determined  by  Equation  2-13. 

When  Equation  6-3  results  were  compared  with  subsonic 
experimental  data,  it  was  deduced  that  better  agreement  would 
be  achieved  if  the  predicted  peak  Cl*  were  to  occur  at  M  <  1.0 
(for  moderate  and  high  AR)  and  if  the  rate  of  increase  in  Cl< 
with  increasing  M  were  larger.  Consequently,  Equation  6-3  was 
altered  to 

(c  \  „  AR/57.3 

Basic 


(6-5) 

where  the  sectional  lift-curve  slope  aD  equals  2 7r,  and  M*  is 
the  limiting  M  for  the  application  of  Equation  6-5.  M*  is  a 
function  of  AR  and _Ac/2  and  is  defined  by 

r  1  2 

M*  -  Mo*  +  (1-Mo*)  [1~(cosAc/2)eJ  (6-6) 

where 

*k*  *  (10  +  0.91AR3)/ (10  +  AR3) 

In  effect,  M*  is  the  Mach  number  at  which  the  rate  of  increase 
in  CLot  with  M  begins  to  decrease.  Note  that,  at  M^O,  Equations 
6-3  and  6-5  are  identical.  However,  it  was  also  found  desirable 
(for  improved  correlation)  to  limit  (coaAc/2)e  to  the  range 
<3.94  cos/tc/2  \  0.  Thus,  in  applying  Equation  6-5,  sweep  angles 
of  less  than  20  degrees  are  treated  as  having  a  value  of  cos^Ac/2  ■= 
0.94. 


f  I 

l+Vl+jl-<cos/l 


<s*> 


n  /  o  «*  ►» 

II 


AR 

2(cos^c/2) 


72 


•*. 


For  thin  wings,  experimental  levels  of  Cl m.  characteristically 
reach  a  peak  at  speeds  somewhat  greater  than  M*.  With  low  sweep 
and  moderate  to  high  AR,  the  peak  occurs  below  M“l,  while  with  high 
sweep  and/or  low  AR,  the  peak  may  be  at  or  above  sonic  speeds.  At 
speeds  well  above  sonic.  Cl*/  then  decreases  with  increasing  M 
and,  when  the  leading-edge  becomes  supersonic  (M  >  cos^le "*),  the 
level  approaches  the  two-dimensional- theory  level  of  Cl^  “ 

4/(57 .3^ M?-"l)  •  To  emulate  these  trends,  Equation  6-5,  was  modi¬ 
fied  by  a  factor  term  and  an  adder  term,  each  to  be  applied  only 
at  M  >  M*.  These  new  terms  are  included  in  the  modified  equa¬ 
tion: 


<eu> 


1/57.3 


Basic 


a/<W(s>  J  ?'/4 


*  2 


M  *>  M* 


(6-7) 


where  Cl,<0  is  defined  as  (CLx)gas^c  at  M  *  M*  in  radians,  and 


w  2 

(S'  -  (M-M*)  [l  +  (M*/M)yJ 
Y  -  (2  -  2/3  -  ^2) 


and 


Z  *  M*C^o  + 


AR 


0  'hAR,  TrAR  1W_  _  , 

3  C?- - i)(cos>{c/2) 

uxo  Wo 


2/3 


Wings  having  thick  airfoils  undergo  a  degradation  in  Cl  <* 
beginning  at  M  >  M*.  The  level  of  versus  M  dips,  usually 

reaching  a  minimum  at  M  <  1.0,  and  then  recovers  to  a  second 
peak  level  at  M  >  1.0.  To  account  for  this  phenomenon,  the  basic 
CL<*  equations  have  been  modified  by  a  factor,  Kt,  as  defined  by 


f4d1(l-  d\)J  3  •  V  ;  MX  /  M  <  M2 

f4(T2(l-<f2)J  3  y  ;  M2  ^  M  f  M3 

Hi  >,  H  M3 


(6-8) 


where  ,  <j~i,  C72>  M^,  M2,  M3  are  as  defined  by  the  equations 
below: 
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predicted  only  if  the  wing  streamwise  airfoil  t/c  er-;  ;jeds  the 
limit  thickness  defined  by  Equation  5-9.  The  limit- thickness 
boundary  was  established  from  the  statistical  boundaries  pre¬ 
sented  in  Reference  24. 

Another  factor  in  the  wing  Cl u  prediction  equation  (Eq.  6-1), 
is  the  fuselage  interference  factor,  defined  by 

Kb  -  (1  +  d/b)(l  -  d/b)f  (6-11) 

where 

16+3AR2 

8+5AR2 

b  =  wing  total  span 

d  -  body  total  width  at  wing  junction 

The  factor  Kb  accounts  for  the  change  in  wing  lift  due  to  the 
body  segment  which  enshrouds  the  wing  and  to  the  wing- induced 
lift  on  that  body  segment.  Based  on  semi-empirical  derivations 
presented  in  Reference  25,  Kb  is  independent  (to  the  first  order) 
of  M.  It  is  noted  that  the  total  lift  of  a  wing/body  configura¬ 
tion  is  derived  by  adding  the  body-alone  lift  to  that  obtained 
for  the  wing-alone  as  modified  by  the  factor  Kb- 

Application  of  Equations  6-2,  6-5,  6-7,  6-8,  and  6-11, 
for  the  prediction  of  (Cl^w-B  yield  trends  as  sketched 

in  Figure  33.  It  is  noted  that  the  technique  is  strongly 
dependent  on  the  value  of  (t/c)  in  the  transonic  speed  regime 
if  (t/c)  exceeds  the  limiting  value  defined  in  Equation  6-9. 

Substantiation  of  the  Ct^..  prediction  technique  described 
above  in  the  form  of  comparisons  with  a  wide  range  of  experi¬ 
mental  data  is  presented  in  References  1  and  26.  The 
derivation  of  Cl for  Mach  numbers  greater  than  M*,  presented  in 
original  Reference  26,  relied  heavily  on  transonic-bump  test  data, 
which  characteristically  produces  a  trend  such  as  shown  in 
Figure  33.  The  new  derivation  (Equation  6-7)  relies  on  sting-mounted 
test  data,  which  produces  a  less  abrupt  transition  in 
(dCL^/dM)  in  the  transonic  region. 


6 . 2  Supercritical  Wing  Lift-Curve  Slope 

The  wing  lift-curve  slope  technique  described  above  in 
Section  6.1  would  underpredict  ClC)<  when  compared  against 
supercritical  wing  data.  In  a  study  of  the  factors  affecting 
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1 


m1+m3 

NOTE:  — = —  is  approximately 


Figure  33  Typical  Lift-Curve  Slopes 
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lift,  it  was  concluded  that  the  reason  Cl  ^  was  being  under¬ 
predicted  is  due  primarily  to  the  thickness  correction 
factor,  Kt  (Equation  6-8),  and  the  incompressible  sectional 
lift-curve  slope,  aQ  (Equation  6-3). 

For  supercritical  wings,  the  onset  of  the  lift  diver¬ 
gence  Mach  number,  Mj_,  is  delayed  to  a  higher  value  as  com¬ 
pared  to  conventional  wings.  For  supercritical  wings,  the 
factors  M]_,  M2,  and  TJ"  in  Equation  6-8  are  modified  as 
follows : 


(M1>SCW  ‘  <mi>conv  +  °-09 

(Mo)  =  (Mo)  +  0.045 
1  SCW  1  CONV 


^sew 


^CONV 


(M3-Ml)SCW 

(M3-M1)conv 


This  modification  delays  the  thickness  correction  factor  to 
a  higher  Mach  number  and  also  decreases  the  extent  by  which 
Kt  is  reduced  at  . 

Supercritical  wings  have  a  higher  sectional  lift-curve 
slope  compared  with  conventional  wings.  The  program  uses 


1  + 


1.174  t/c 

V7^2 


for  supercritical  wings  in  place  of  a0/27r-  1  for  conventional 
vin^s • 


6.3  Tail  Lift-Curve  Slope 

The  lift-curve  slope  of  the  horizontal  tail  can  be 
estimated  from  the  equation 

(CL*>t  *  <CL*  >Tj%(B)+KB(W)]  (1-  J£>  s~  (6~13> 

where  (CjJ^-j,  is  the  exposed-area  lift-curve  estimate  for  the 
tail;  and  Kg^  are  the  Pitts,  Nelson,  and  Kaattari  body- 

lift  carry-over  factors  (Reference  27);  is  the  downv7ash 

gradient;  qt/q^  is  the  dynamic  pressure  ratio;  and  Sj  is  the 
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exposed  tail  area.  The  exposed-area  lift-curve  for  the  tail  is 
estimated  by  use  of  the  exposed  planform  of  the  tail  and  the 
method  described  in  Section  6.1. 


t 

\ 


6.3.1  Downwash  at  the  Tail 


An  empirical  method  of  estimating  the  low-speed  downwash 
gradient  behind  straight- tapered  wings  is  given  in  the  DATCOM  by 


4.44 


1.19 


(6-14) 


The  factors  KA,  and  Kjj  are  wing  aspect  ratio,  wing  taper 

ratio,  and  horizontal- tail- location  factors,  respectively, 
determined  from 


Ka  =  1/AR  -  l/Cl+AR1'7) 


K  A 


10-3A 

7 


and 


K, 


H 


d-Mb) 

(2yt'/b)1/3 


where  AR  and  are  the  wing  aspect  and  taper  ratios,  respectively, 
ht  is  the  height  of  the  tail  relative  to  the  wing  chord  plane  and 
/t’  is  the  distance  between  the  exposed  MAC  of  the  wing  and  the 
exposed  MAC  of  the  tail.  At  higher  speeds  the  effect  of  compres¬ 
sibility  on  downwash  is  approximated  by 


36  =  /-*6\  (cLrf)M 

i^)0  (cOo 


(6-15) 


where  ( Cl^) 0  and  (cL<x)m  are  t^ie  lift:- curve  slopes  at  low 

speed  (M=0. 1)  and  at  the  appropriate  Mach  number,  respectively. 
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6.3.2  Dynamic  Pressure  at  the  Tall 


The  method  for  estimating  the  dynamic-pressure  qt/q  .*  at 
ihe  tail  is  based  on  the  DATCOM  method  which  relates  the 
d  namic-pressure  ratio  to  the  drag  coefficient  of  the  wing. 

The  steps  involved  in  determining  the  dynamic  pressure  at  some 
distance  aft  of  the  wing  root  chord,  outlined  in  Section  4.4.1 
of  the  DATCOM  report,  are  as  follows 

1.  Calculate  the  half-width  of  the  wing  wake  by 


~  =  0 • 68  v Cdq  (£  +  0.15) 


(6-16) 


where  x  is  the  longitudinal  distance  measured 
from  the  wing-root-chord  trailing  edge,  Zw  is 
the  half -width  of  the  wake  at  any  position  x, 
and  Cdo  is  the  wing  zero-lift  drag  coefficient. 

2.  Calculate  the  downwash  i  .i  the  plane  of  symmetry 
at  the  vortex  sheet  by 


i.bZ 

7MR 


(CL«2) 


3.  Determine  the  vertical  distance  Z  from  the 
vortex  sheet  to  the  quarter-chord  point  of 
the  MAC  of  the  horizontal  tail  by 

Z  =  x  tan  ( y +  6  -  2) 

where  7  =  tan-1-  (ht/^t). 


(6-17) 


(6-18) 


4.  Determine  the  dynamic-pressure-loss  ratio  at 
the  wake  center  by 

/2)q.  _  2.42(CDo)1/2 


(-“>  = 
v  q  Jo 


(£  +  0.3) 


(6-19) 


5.  Determine  the  dynamic-pressure-loss  ratio  for 
points  not  on  the  wake  centerline  by 
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(6-20) 


4a 


2  ,7T 

COS  (tj- 


2  Z 


w 


6.  Determine  the  dynamic  pressure  ratio  at  an 
arbitrary  distance  x  aft  of  the  wing-root- 
chord  trailing  edge  by 


at  =  L  _  A<? 


(6-21) 


6 , 4  Body  Lift-Curve  Slope 


As  shown  in  Reference  28,  the  linearized  lift-curve  slope 
for  a  body  can  be  expressed  as 


_F _ 

sRef 


(6-22) 


where  L  is  the  body  length,  Lfl  is  the  forebody  length,  F  is  the 
body  cross-sectional  area,  and  k^  is  a  linear-potential  lift- 
curve-slope  parameter.  The  factor  k^  (a  function  of  body  width 
b,  body  height,  h,  and  the  perimeter  of  an  elliptical  body  with 
equal  area,  p)  is  determined  from  the  curve  given  in  Figure  34, 

Uiuui  Is  L<jK.eu  JLJ.UUI  n.ejLCicii(.c 


6.5  Angle  of  Attack  at  Zero  Lift 

The  angle  of  attack  at  zero  lift,  aLo>  *-s  determined  from 

aLo  ~  ^“Lo^cAMBER  +  ^^TWIST  +  INCIDENCE  (6_23) 

The  effect  of  camber,  Cl^  on  ot^  is  calculated  from 

(6-24) 


.  a«Lo 

aLo)  CAMBER  ^dCL  Ld 


is  shown  plotted  in  Figure  35,  which  is  obtained 
onal  data. 


where  (dofLo/SC^  ,)  i 
from  two -dimens  Iona 

The  increment  in  due  to  wing  twist,  r  ,  is  calculated 


from 


/  v  _  /  d«Lo,  _ 

^Lo^  TWIST  ^ 


(6-25) 
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where  5C*Ln/a  T  “  °'093  "  0.00^571^+  0.5761^-  0.2645^  and 

tan-1  (tarU^/4/^)  (deg)  The  equation  for  (O^Lo^'t  )  waS 
obtained  from  a  curve  fit  of  the  parametric  data  reported  by 
Gilman  and  Burdges  (Reference  29)  for  wings  with  linear-element 
twist. 

The  angle  of  attack  in  the  program  is  measured  relative 
to  the  wing  root-chord  reference  plane.  For  variable-sweep 
configurations  the  angle  of  attack  for  any  sweep  position  is 
measured  from  the  wing  chord  plane  in  the  forward  sweep  posi¬ 
tion.  The  increment  in  O^lq  due  to  wing  and  horizontal  tail 
incidence  is  calculated  from 

)  <iw)  +  (CU)  Uw-it> 

(c(  V  =  BODY _ TAIL _ 

1x5  INCIDENCE  CL* 


+  (Vt 


WREFy 


(6-26) 


where 


(CI^) 


Rnnv 


(CL*> 


MnCT7 


~W(B) '~B(W)  W'J 


When  M  >  1,  the  contribution  of  camber  and  twist  to 
is  set  equal  to  zero  and  only  the  incidence  effect  is  continued 
supersonically. 


Nonlinear  Lift  of  High- Aspect-Ratio  Wii 


The  lift  characteristics 
of  a  high-aspect-ratio  wing 
is  illustrated  in  the  sketch. 

A  high  aspect  ratio  is  defined 
as  AR)  ARlqw>  where  ARj/jy  is 
defined  in  Section  6.9  (Equa¬ 
tion  6-37).  The  lift  varies 
linearly  with  angle  of  attack 
up  to  Cls>  after  which  the 
lift  variation  becomes  non-linear. 
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The  angle  of  attack  for  a  specified  lift  coefficient  is 
calculated  from 


where 


CL 


+  +AC/- 


Cl-  '  Lo 


0; 


(6-27) 


CL^  CLS  or  M  >  l-° 
2 


A  (X 


/  CL-CLs  y 

(CLmax  CLs)  max;  ^  ^  Cl  ^  CL|nax 

^  (CLmax-CLs)  '  ^max  +  5  ?  °L  >  CW 


(6-28) 


and  CLS  “  Lo  '  2A^max) 


The  prediction  of  Cl^x  and  g^p(max  f°r  high-aspect-ratio  wings 
is  discussed  in  Section  6,9. 


6.7  Nonlinear  Lift  of  Low-Aspect-Ratio  Wings 


The  subsonic  charac¬ 
teristic  of  a  low-aspect- 
ratio  wing  is  illustrated 
in  the  sketch.  The  total 
lift  is  equal  to  the  po¬ 
tential  lift  plus  the  vor¬ 
tex-induced  lift  from  the 
leading  edge  and  tip  of 
the  wing.  The  equation 


for  lift  can  be  expressed  as 


CL  *  CLP  +  CLV  +  CLTv  (6-29) 

The  method  of  predicting  each  of  these  terms  is  discussed  in 
the  following  subsections. 


6.7.1  Potential- Flow  Lift 


The  potential -flow  lift  it,  determined  from 

Clp  =  Kp  sin  C<cos^o(d-  (6-30) 

where  Kp  is  the  lift-curve  slope  given  by  small-disturbance 
potential- flow  lifting-surface  theory,  and  the  trigonometric 
terms  account  for  the  leading-edge  separation  effects  (Refer¬ 
ence  30)  .  The  value  of  Kp  is  the  lift- curve  slope  (Cl«),  con¬ 
verted  to  radians;,  obtained  from  Equation  6-1.  The  factor  Clq 
is  the  lift  at  zero  angle  of  attack  predicted  by 


CL0  ~  -Cl*  0(1^ 


6.7.2  Leading-Edge  Vortex  Lift 

The  leading-edge  vortex  lift  is  determined  from 

Cl^  '  (1-R)  •  FVB  '  Ky  sin^oCcosoC  (6-31) 

In  this  equation,  developed  in  Reference  31,  the  sharp- leading- 
edge  suction  analogy  of  Polhamus  is  modified  to  account  for 
round-leading-edge  and  vortex -breakdown  effects. 

In  Figure  36  (taken  from  Reference  36)  the  theoretical 
variation  of  the  vortex  lift  factor  Ky  with  aspect  ratio  and 
cutout  factor  is  shown.  The  factor  Fyg  (shown  in  Figure  37) 
is  a  vortex-breakdown  factor,  which  was  obtained  from  the 
ratio  of  experimental  data  to  theoretical  for  sharp- leading- 
edge  delta  wings.  The  factor  R  in  Equation  6-31  is  a  lead¬ 
ing-edge  suction  parameter  (Reference  33).  For  a  sharp 
leading  edge,  the  suction  parameter  is  near  zero;  for  a 
rounded  leading  edge,  the  suction  parameter  is  near  unity 
at  low  alphas.  The  variation  of  R  versus  oL  is  shown  in 
Figure  38  as  a  function  of  thickness  ratio. 
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6.7.3  Tip-Vortex  Lift 


\ 


For  wings  having  a  tip  chord  greater  than  zero,  a  tip  vortex 
forms  that  induces  an  additional  lift  contribution  on  the  wing. 

At  low  angles  of  attack  the  flow  around  the  wing  leading  edge  and 
tip  is  attached,  and  a  vortex  sheet  is  formed  at  the  trailing  edge 
(Figure  39a).  At  slightly  higher  angles  of  attack  f figure  39b), 
the  flow  possibly  will  make  the  turn  around  the  leading  edge  of 
the  wing  without  separating,  but  the  flow  around  the  tip  separates. 
In  this  stage,  the  flow  forms  a  vortex  sheet  consisting  of  a  hori¬ 
zontal  part  originating  from  the  trailing  edge  and  two  vertical 
sheets  attached  to  it  originating  from  the  two  sides  of  the  wing. 
Kuchemann  (Reference  34)  noted  that  a  spanwise  cross-section 
through  the  vortex  sheet  has  the  same  shape  as  that  obtained 
behind  a  wing  with  end  plates.  The  height  of  the  "end- plate 
vortex"  or  tip  vortex  is  approximated  by 


h/b 


Of  ^T  JL 
2  c  AR 


(6-32) 


where  Op  is  the  tip  chord.  With  the  height  of  the  tip  vortex 
known,  the  incremental  tip-vortex  lift  can  be  expressed  as 


Z+  n/T+Z2 
xZ+  /  l-f-(xZ)  ^ 


Krv  sin  a 

t' 


(6-33) 


where 


Z  -  2  cosAc/2/AR 

x  =  1. 0014-1. 969(h/b)+3.0021(h/b)2-2.0072(h/b)3  (6-34) 

Equation  6-33  was  derived  (Reference  34)  by  modifying  the  Helmbold  : 

lift  equation,  where  the  effect  of  the  end-plates  are  expressed  i 

as  a  factor  1/x  to  the  aspect  ratio.  Equation  6-34  is  a  curve  ! 

fit  of  the  end-plate  effect  shown  in  Figure  III. A. 4-1  of  Refer-  \ 

ence  4.  | 

It  can  be  seen  from  Equation  6-32  that  the  end-plate  effect 
becomes  smaller  with  increasing  aspect  ratio.  This  explains 
why  the  end-plate  effect  of  the  tip  vortex  has  rarely  been 
noticed  for  wings  of  moderate  and  large  aspect  ratios, 
although  it  always  existed.  The  tip-vortex  method  is  shown 
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in  Reference  34  to  give  good  agreement  with  experimental  data 
for  unswept  rectangular  wings  ranging  in  aspect  ratio  from  0.5 
to  2.0.  It  is  stated  in  Reference  34  that  the  end-plate 
analogy  can  be  used  for  straight  or  swept  wings . 

The  nonlinear  lift  calculated  by  Equation  6-29  is  limit  ed 
by  the  maximum  lift  coefficient,  Ci_aX,  of  the  wing.  The  value 
of  and  the  angle  at  which  maximum  lift  occurs, oCmaxi  are 

predicted  by  the  low-aspect-ratio  method  given  in  Section  6.9. 
The  leading-edge  vortex  is  then  limited  by  the  condition 

ClV  '  ClMAX  “  (ClP*max  "  <CLtv)MAX  (6-35) 

where  (Clp)maX  anc*  (^Ltv^MAX  are  t^e  values  calculated  for  the 
potential  lift  and  the  tip-vortex  lift  at  the  maximum- lift  angle 
of  attack.  If  (Cl^O^^X  ^  (Cl)max>  it  is  assumed  that  the 
leading-edge  and  tip  vortices  are  too  weak  to  add  much  lift  to 
the  wing,  and  the  high-aspect-ratio  method  discussed  in  Section 
6,6  is  then  used  to  predict  the  lift  up  to  stall. 

Results  of  applying  the  nonlinear  lift  prediction,  procedure 
are  shown  in  Figures  40  through  43.  The  data  were  taken  from 
Reference  35,  which  reports  on  a  test  of  a  series  of  clipped 
delta  wings.  The  program  results  are  shown  as  the  solid  lines 
for  the  complete  lift  and  as  dashed  lines  for  the  initial  value 
of  the  lift-curve  slope,  Cl,*..  The  results,  in  general,  are 
good  and  indicate  a  substantial  amount  of  the  lift  is  due  to  the 
vortices . 


6 . 8  Nonlinear  Lift  of  Cranked  Wings 

The  method  available  in  the  Large  Aircraft  program  for 
predicting  the  subsonic  lift  variation  of  cranked  wings  is 
based  on  the  technique  presented  in  Reference  9.  This  method 
assumes  that  as  the  outboard  panel  of  a  cranked  wing  experiences 
stall,  the  inboard  panel  still  continues  to  lift.  This 
behavior  is  believed  to  be  caused  by  the  influence  of  the 
leading-edge  vortex  of  the  inboard  panel.  Consequently,  the 
flow  field  is  similar  to  that  of  a  low-aspect-ratio  delta 
wing  with  leading-edge  separation. 

The  method  employs  the  results  of  a  data  correlation 
accomplished  to  provide  a  technique  for  determining  the 
nonlinear  lift  of  double-delta  wings.  It  is  hypothesized 
that  the  nonlinear  lift  of  a  cranked  wing  should  be  similar 
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to  that  of  a  double-delta  wing.  The  nonlinear- lift  curve 
construction  technique  for  cranked  wings  is  shown  in  Figure 
44. 


Figure  44  Construction  of  Nonlinear  Lift  Curve  for 
Cranked  Wings 


Except  for  some  slight  refairing  of  the  data  correlation 
curves  to  account  for  aspect  ratios  of  less  than  1.0,  the  method 
employed  in  the  program  is  essentially  the  same  as  that  presented 
in  Reference  9.  The  nonlinear  angle  of  attack  for  a  given  Cl  is 
calculated  from 


CX  = 


I  Ai 
a  Cl*?0 


l/n 


(6-36) 


where 


Cl  is  the  linear  lift-curve  slope  (including 
outboard  panel) 

Ai  is  the  aspect  ratio  of  the  inboard  panel, 
defined  as  shown  in  Figure  45. 
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is  the  nondimensional  spanwise  ordinate 
for  the  break  point  in  the  cranked  wing 

a,n  are  correlation  constants  (shown  in 
Figure  46}  derived  by  modifying  the 
Reference  9  method. 

In  the  program  it  is  assumed  that  o(  stall  is  angle  of 

attack  corresponding  to  CijDg,  as  shown  in  Figure  13  of  Section  4. 
The  accuracy  of  the  nonlinear  angle-of-attack  prediction  tech¬ 
nique  is,  of  course,  strongly  dependent  on  knowing <X stall* 


6.9  Maximum-Lift  Coefficient 


The  method  used  in  the  Large  Aircraft  program  to  estimate 
the  maximum  lift  and  angle  of  attack  for  maximum  lift  is  based 
on  the  DATCOM  method  for  both  the  low-  and  high-aspect-ratio 
wings  (Reference  2).  The  maximum  lift  of  high-aspect-ratio 
wings  at  subsonic  speeds  is  directly  related  to  the  maximum 
lift  of  the  wing  airfoil  sections.  The  wing  planform  shape  is 
a  secondary  influence  on  the  maximum  lift  obtainable.  However, 
for  low-aspect-ratio  wings,  the  wing  planform  is  the  primary 
effect  on  maximum  lift, while  sectional  characteristics  are 
secondary.  The  program  uses  the  criteria  established  in  the 
DATCOM  method  by  the  equation 


AR 


L0W  (Ci+1)  cosAle 


(6-37) 


where  C^  is  a  function  of  t-aper  ratio,  as  given  in  Figure  54. 
If  AR  >  ARlow,  the  high-aspect-ratio  method  is  used,  and 
if  AR  ^  ARlOW»  the  low-aspect-ratio  method  is  used.  These  two 
methods  are  described  in  the  following  subsections. 


6.9.1  High-Aspect-Ratio  Method 

The  DATCOM  method  is  an  empirically  derived  method  based 
on  experimental  correlations  of  high-aspect-ratio,  untwisted, 
and  constant-section  wings.  The  equations  for  maximum  lift 
and  the  angle  of  maximum  lift  are  as  follows: 


cLmax  ( 


C4iAX 

c/max 


^  c/max  +^ClMAX 


(6-38) 
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(6-39) 


OC  MAX 


rCLMAX 
CL 


^^NAX 


The  first  term  in  Equation  6-38  is  the  maximum  lift  coefficient 
at  M-0.2;  the  second  term  is  the  lift  increment  to  be  added  for 
Mach  numbers  greater  than  0.2. 

The  factor  (Cl^ax/C/j^x)  is  computed  by  a  curve  fit  of  the 
curves  in  Figure  4. 1.3. 4- 14a  in  the  DATCOM  given  by 


lMAX 

^MAX 


-  A  -  B  Ay' 


(6-40) 


where 


0; 


Ay'  =  J^y-i.4; 


1.1: 


Ay  c  1.4 
1-4  <  2.5 

y  >  2.5 


and  the  terras^  A  and  B  are  plotted  in  Figure  47  as  a  function 
of  sweep.  (The  Ay  is  defined  by  Equation  2-16.) 

The  increment  to  Clj^x  due  bo  Mach  number  is  computed  from 
a  curve  fit  of  the  curves  of  Figure  4.1.3.4-15  in  the  DATCOM 
given  by  ^ 


ci-max  - c  +  <D-C><sr> 


(6-41) 


where  the  terms  C  and  D  are  plotted  in  Figure  48  as  a  function 
of  Ay  and  Mach  number. 


The  section  maximum  lift  coefficient  at  M*0.2,  C/may  is 
computed  from 


C^MAX  ^iMAX^  +  ^%1AX  (6-42) 

Base 

where  (C>^^x)gase  and  C/mAX  are  shown  plotted  in  Figures 
49  and  50  as  a  function  of  the  sharpness  factor,  maximum- thick¬ 
ness  location,  and  camber. 
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The  angle -of -attack  increment  for  maximum  lift, 
is  obtained  from  Figure  51,  which  is  taken  from  Section  4. 1.3.4 
in  the  DATCOM. 


6.9.2  Low-Aspect-Ratio  Method 

The  empirical  equations  in  DATCOM  for  estimating  subsonic 
maximum  lift  and  angle  of  attack  for  untwisted  low-aspect-ratio 
wings  are 

ClMAX  *  ^MAX^ase  +  *^ClMAX  (6-43) 

4XMAX  =  (<*MAx)Base  +  ^MAX  (6-44) 

The  base  value  of  CLj^yj  is  obtained  from  Figure  52  if  the 
position  of  maximum  airfoil  thickness,  Xx»  is  forward  of  the 
35-percent  chord  point,  and  from  Figure  53  if  Xj-  is  aft  of 
the  35-percent  chord.  The  values  of  C]_,  and  C2  are 

obtained  from  Figure  54,  the  base  Qj^ax  from  Figure  55,  and 
the  value  of^iWMAX  from  Figure  56.  (Figures  52  through  56 
are  taken  from  Section  4. 1.3.4  of  the  DATCOM.) 


6.9.3  Tail-Lift  Contribution  to  Cl>[ax 

Because  the  horizontal  tail  usually  has  a  smaller  aspect 
ratio  compared  to  the  wing  and  is  In  a  downwash  field  that 
counteracts  the  effect  of  angle  of  attack  to  some  extent,  it 
is  assumed  that  the  tail  does  not  stall  before  the  wing.  The 
lift  generated  by  the  tail  at  the  angle  of  attack  of  wing  stall 
is  added  to  the  wing  maximum  lift  coefficient  to  obtain  the 
configuration  maximum  lift  coefficient.  The  configuration 
maximum  lift  is  given  by 

CLmax  ‘  (CI-MAX)U.  +  WW.  ..  <6-45 

Wing  Tail 


where 


^"'Lmax'  ^  ’  57.3  sin  ^max(cos  0(max^ 

The  term  (Cl^)^  is  the  lift-curve  slope  of  the  horizontal  tail 
as  determined  in  Section  6.3. 
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Figure  52  Subsonic  Maximum  Lift  of  Low-Aspect-Ratio  Wings  (  X^,  <  35%) 
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Figure  53  Subsonic  Maximum  Lift  of  Low- Aspect-Ratio  Wings  (xt  5  35%) 
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7 .  MDMENT 


The  moment  of  a  wing- body- tail  configuration  can  be  repre¬ 
sented  by  the  equation 

,  Xac  cRe 

Cui  “  Cmo  +  ^n"  Cr^  T~  ClWB  ■  CLtAIL  *  ^Ht/c  (7_1) 

where  n  is  the  chordwise  distance  to  the  moment  reference  point 
measured  in  exposed  wing  root  chord  (Cr^,  Xac/CRe  is  the 
aerodynamic  center  location  relative  to  Cr6,  and  the  last  term 
represents  the  moment  contribution  of  the  tail  lift  times  the 
moment  ar®,{HT,  determined  from  Equation  (2-9). 

The  following  subsections  discuss  the  method  used  to  predict 
the  elements  in  Equation  7-1  along  with  a  method  of  determining 
trimmed  lift  curves  and  drag  polars. 


7.1  Zero-Lift  Moment 


The  method  of  predicting  the  zero-lift  pitching  moment  for 
a  wing-body  configuration  considers  only  the  effect  of  the  wing 
°n  Cipq  and  does  not  include  the  effect  of  an  assymmetrical 
fnse!;=;ge  or  the  effect  of  stores  and  nacelles  located  near  the 
Wi'g.  However,  the  Cmo  prediction  method  in  the  Large  Aircraft 
program  can  be  adjusted  by  input  so  as  to  match  the  test  data 
Cmo  on  a  similar  configuration. 

The  subsonic  zero-lift  pitching  moment  for  wings  with 
linear  twist  ,  up  to  the  critical  Mach  number,  is  given  in  the 

A  mni  .  w  _  _  *  CD 


T%  A  rPr*i  t\A  -  ~ 
i/niwn  ao 


=1' 


(S>>r.0  ♦  r  ! 

Z  >  Vl  -  M^OS 


(7-2) 
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where  (Cmo)^=:g  is  the  Cn^  of  an  untwisted  wing  and  (AQno/r)  is 
the  change  in  wing  zero- lift  pitching  moment  due  to  a  unit 
change  in  wing  twist,  'f.  The  parameter  (/iC^ It)  was  obtained 
from  lifting- line  theory  and  is  shown  in  Figure  57. 


The  Cn^  of  an  untwisted  wing  is  obtained  from 

ARcos3tc/4 

(CmoV=0  “  AR+2 co sAc/4  (C^o)SECT 


(7-3) 
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where  (Cmo)^^  is  the  average  section  pitching  moment  coef¬ 
ficient  determined  by  averaging  the  section  Cmo  for  each  wing 
panel,  using 


(Cmo) sect 


m 


N 

Z/s  t 

i=l 


^/(c/d>iSi 


where  (Qn/Cy^)  is  the  theoretical  pitching  moment  divided  by 
the  section  design  lift  of  the  airfoil  camber  line  obtained 
from  Table  4.1.1-D  in  the  DATCOM. 


7 . 2  Aerodynamic  Center 

The  aerodynamic  center  location  of  a  wing-body  configu¬ 
ration  is  given  in  the  DATCOM  as 


Aac 

CR. 


(C°w(b>*(^Jb(w)  (ClJb(w> 


(  CL  «.)  fi+'  C“*r  W  (B)+(  )  B  (  W) 


(7-4) 


where  the  Xac/CRe  terms  are  the  chordwise  distances  measured 
in  exposed  wing -root  chords  from  the  apex  of  the  exposed 
wing  to  the  aerodynamic  center,  positive  aft.  The  sub¬ 
scripts  N,  W(B),  and  B(W)  refer  to  the  lift  and  aerodynamic 
center  contribution  of  the  forebody,  exposed  wing,  and  the 
wing- lift  carryover  on  the  body,  respectively. 


7.2.1  Aerodynamic  Center  of  Forebody 

The  subsonic  location  of  the  aerodynamic  center  for 
forebodies  with  ogive  nosecones  is  approximated  in  the 
DATCOM  as 

Ljj+l.eCXLg.Ljj) 

CRe  (7 

Figure 58  defines  the  geometric  parameters,  Ln»  Xle,  and  CRe. 
The  supersonic  forebody  aerodynamic  center  is  obtained  from 
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(7-6) 


/Xac\  =  ^LE  / XCP  _  1\ 

\CRe/N  CRe  \  Jb  '  / 

where  the  term  X^p/^g  is  obtained  from  Figure  59  (Figure  4. 2. 2.1- 
23a  in  DATCOM). 


7.2.2  Aerodynamic  Center  of  Wing  (Trapezoidal,  Single  Panel) 

The  aerodynamic  center  of  the  exposed  wing  is  determined 
from  the  DATCOM  charts  presented  in  Figures 60a  through  60f. 
These  charts  are  valid  for  subsonic  Mach  numbers  less  than  Mach 
critical  and  supersonic  Mach  numbers  greater  than  1.2.  For 
transonic  conditions,  the  data  presented  in  the  DATCOM  in  terms 
of  transonic  similarity  parameters  (Figures  61a  through61d)  are 
used  to  determine  the  aerodynamic  center  position. 

The  procedure  for  obtaining  aerodynamic  center  can  be 
summarized  as: 


For  M  <  McR; 


I 


For  Mcr+,05  M  ?  ^CR, 


xac  \  /  xac  \  /xac  \ 

:crJ~{crJ  CRJ  -  [CrJ 


M-Mcr 

.05 


For  \/  l+(t/c)2/3  >  M  >Mcr  +  .05, 


(7-7) 


(7-8) 


(7-9) 


(7-10) 
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Ref.  2,  FIGURE  4.2.2. 1-23  a 
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Ref.  2,  FIGURE  4.1. 4.2 -22 
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Figure  60  Wing  Aerodynamic -Center  Position 
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(7-11) 


For  M  M.2, 


where  (Xsc/CRe)'  is  read  from  Figure  60  and  (Xac/CRe)"  is  read 
from  Figure  61. 


7.2.3  Aerodynamic  Center  of  Wing  (Cranked  or 
Double  Delta) 

The  prediction  of  the  a.c.  location  of  cranked  or  double¬ 
delta  wings  is  taken  from  the  method  developed  at  Convair  Aero¬ 
space  as  reported  in  Reference  10.  The  non-straight-tapered 
wing  is  divided  into  two  panels,  with  each  panel  having  conven¬ 
tional,  straight- tapered  geometry.  The  individual  lift -curve 
slope  and  a.c.  are  estimated  for  each  panel,  using  the  technique 
described  above  for  the  trapezoidal  wing  and  treating  each  con¬ 
structed  panel  as  an  isolated  wing.  The  individual  lift  and 
a.c.  location  for  each  constructed  panel  are  then  combined,  using 
an  inboard-outboard  weighted-area  relationship 


^ac 

C*e 


(Cl^Sl  +  <.CLj0S0 


(7-12) 


where  the  outboard  wing  a.c.  is  referenced  to  the  inboard  root 
chord  length  given  by 


/Xa„\  /X3C\  Cr 

•  ■  I  ^  V  I  U 

V <W0  =  l %/0 


(tanALE>o 


th/9\  . 


The  geometry  for  the  inboard-outboard  panel  arrangement  is 
illustrated  in  Figure  62. 


7.2,4  Aerodynamic  Center  of  Wing-Lift 
Carryover  on  Body 

The  location  of  the  a.c.  due  to  the  wing- lift  carryover  on  the 
body  is  determined  by  use  of  the  DATCOM  method.  For  /?ARe  ^  4  the 
subsonic  a.c.  location  is  obtained  from 


1  ,  fc-d 
5  +  2CT 


tan  A 


f(d/b) 


(7-13) 


B(W) 


where  ARe  is  the  expos ed-wing  aspect  ratio  and  the  factor  f(d/b)  is 
shown  plotted  in  Figure  63.  For  y9ARe  <  4  the  a.c.  location  is 
determined  from 


(7-14) 


kij&imm  +ft)’ 


where  (X^c/Cr  )’  i-s  the  a.c.  location  determined  from  Equation  7-13 
and  (Xac/'CRe)R  is  the  theoretical  location  for  /?ARe=0  determined 
from  the  equation 


/Xac\"  _  1 
\CRe  /  "  8 


ARe(l+Xp)  tan  A 


(7-15) 


where  Ae  is  the  exposed-wiug  taper  ratio.  Equation  7-15  is  limited 
to  values  less  than  or  equal  to  0.5.  For  supersonic  conditions  the 
a.c.  location  is  estimated  from  Figure  64,  For  transonic  conditions 
the  a.c.  location  is  determined  by  linear  interpolation  of  the  a.c. 
values  determined  at  the  critical  Mach  number  and  Mach  1.1. 

For  complex  wing  planforms  the  equivalent  wing-sweep  values  are 
used  in  the  subsonic  and  supersonic  a.c.  location  methods. 


7 . 3  Effect  of  Trim  Deflection 

The  effect  of  trim  deflection  can  be  estimated  by  predicting  the 
incremental  change  in  lift,  drag,  and  moment  due  to  tail  deflection 
along  constant  angles  of  attack.  The  total  wing-body- tail  lift, 
drag,  and  moment  can  be  represented  by 


CL  "  CLWB  +  cLat-(«-  «ot)  +  CL*  * 


(7-16) 
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Ref.  2,  FIGURE  43-2.1-35 


rigute  uj  rarameter  Used  in  Accounting  for  Wing-Lift 
Carryover  on  the  Body 


Ref.  2,  FIGURE  4.3.2.1-37 


EFFECTIVE  RATIO  OF  BODY  DIAMETER  TO  ROOT  CHORD,  — 


Figure  64  Aerodynamic -Center  Locations  for  Lift  Carry¬ 
over  onto  Body  at  Supersonic  Speeds 


(7-17) 


CD  *  CDmIN  +  <CI>L)  +  (CDL> 

TOTAL  ^  HT 


+  -  |cLat(a  -Cf0t)  +  CLa  6  ht  4t/C  <7- 

L  W-B  ' 

where  the  induced  drag  of  the  tail  is  predicted  by  use  of  a  tail- 
induced  drag  factor  times  the  square  of  the  lift  generated  by 
the  tail,  i.e. 


18) 


(cDl)ht  -  Kt  j^CLflt(a-a0t)  +  CLa  3HtJ 


(7-19) 


From  Equations  7-16  and  7-17,  the  incremental  change  in  lift  and 
drag  at  constant  a  due  to  a  trim  deflection  can  be  determined  by 


ACL  =  Cr 


5 


HT 


(7-20) 


ACD  =  Kt[cLat(«-<^t)4CL<*  aHT|2  -  Kt  [cL<<t(«-aot)J2  (7-21) 


Equation  7-21  can  be  reduced  to 

ACj)  “  S^HT  +  b^HT^a  ~  ^ot^ 
where  a  =  KtCL  2 

cL 


(7-22) 


from 


b  =  2KtCLat  ^ 

The  factor  Cl 4  for  an  all-movable  horizontal  tail  is  predicted 


CL*  “  (cLa)t  %(B)  ’  ^t'^, 


(7-23) 


where  (0^),.  is  the  exposed-area  lift-curve  estimate  for  the  tail,  and 
^W(B)  1®  tbG  surface  lift  in  the  presence  of  the  body  factor.  The 
inauced  drag  factor  for  the  tail  is  determined  by  the  leading-edge 
suction  method  described  in  Section  4,  where 


Kr  = 


1-Rt  Rt  Sref 

+ 


t  57  •  3KW(  b/cL^)l-  7TARht  (Sex) 


(7-24) 


HT 
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The  Large  Aircraft  program  can  predict  the  lift,  moment, 
and  drag  for  a  fixed  tail  setting  or  the  program  can  solve  for 
h HT  to  trim  out  the  moment .  The  tail  deflection  required  for 
trim  is  obtained  by  setting  Equation  7-18  equal  to  zero  and 
solving  for  £>HT* 


8.  HIGH™ LIFT  SYSTEM  AERODYNAMICS 


The  empirical  methods  for  predicting  lift,  drag,  and 
moment  of  an  airplane  with  flap  and  slats  deployed  consist 
of  adding  the  incremental  effects  of  the  high- lift  system 
to  the  clean  airplane  aerodynamics.  Figure  65  illustrates 
the  manner  in  which  the  incremental  effects  of  a  flap  can 
be  applied  to  the  clean-wing  aerodynamics.  The  following 
techniques  for  estimating  these  increments  were  derived 
from  DATCOM  and  References  36  and  37. 


8. 1  Lift  of  High-Lift  Devices 

The  untrimmed  equation  for  lift  can  be  expressed  as 


CL  -  %  +  cL^(*  -A«) 

The  term  Abi  varies  with  CM  according  to 


(8-1) 


<*-  ^<<*max  -  2*0(max) 


Ao(  = 


foL  "<0^n,!»v+2Ao(  v\  2 


(8-2) 


(c*r 


max'^&max) 


where 


c/ 


max 


CWClo 


+  ^  O^max 


The  clean-wingAc/max  determined  in  Section  6.9.1  is  also  used 


for  the  case  with  flaps.  The  increment  in  Clo  and  caused 

K.r  O  1  - PI  /Icf  1  n  on  Vo  o 

L/  J  u  A.  JLU^/  UWP  WUik  k/W  W  W  k.  ■  >  IU  kV  VI  k«.  V  111 

(8-3) 

AcW  -  AC/max  KCK„  cos7tc/2  (8-4) 


where  Kc  and  are  shown  plotted  in  Figures  66  and  67.  The 
factor  Cl*/ Clot.  ,  determined  from  the  Polhamus  lift  equation 


(8-5) 
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converts  the  sectional  values  ACjq  and  .AC/max  to  the  three- 
dimensional  case.  The  sectional  values  are  either  obtained 
from  input  to  the  program  or  generated  internally  for  certain 
types  of  high- lift  system.  The  method  of  generating  the  section 
values  is  discussed  in  Section  8.4. 

The  increment  in  lift  at  zero  angle  of  attack  is  approxi¬ 
mately  zero  when  a  slat  is  deflected.  The  slat  acts  to  delay 
separation  from  the  wing  leading  edge  and  thus  allows  higher 
angles  of  attack  and,  consequently,  higher  values  of  maximum 
lift  before  the  wing  stalls.  An  estimate  of  the  increase  in 
maximum  lift  of  a  slat  is  represented  by 


where  the  partial  span  effect,  Kg»  is  shown  plotted  in  Figure  68 


8 . 2  Drag  of  High- Li ft  Devices 

The  un trimmed  equation  for  drag  can  be  expressed  by 

CD  -  cD„,in  +K(CL-i4CL)2  (8-7) 

Total 


The  drag-due- to- lift  factor,  K,  with  high- lift  devices  car.  be 
estimated  by 


K  ^Clean 


sRef 

fs-rxT  _ 

^  rj.au 


where  Kciean  tne  clean-airplane  polar  factor  and  Sp^an  is 
the  new  planform  area  of  the  wing  if  the  flap  extends  the 
chord  of  the  wing.  The  drag  due  to  lift  factor  remains 
unchanged  until  a  >  amax  -  2  •  Aamar  after  which  K 
approaches  the  zero-suction  value  given  by  Cdl  =  CL  tan<^max* 

Total  minimum  drag  for  the  high- lift  configuration  is 
expressed  as 


c3min  =  CDmin  +  C»L.G.+  Cl>Flaps+  CDStat.>+C'>I  (8'8) 
Total  Clean 


where  increments  are  summed  for  the  landing  gear  and  the  profile 
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drag  caused  by  the  flaps  and  slats.  The  increment  in  profile 
drag  of  the  flaps  and  slats  can  be  estimated  from  sectional 
drag  data,  using 

^CDpiap  =-ACdf  cos71hl  *  Kd  (8"9) 


^cDSLat  =^Cds  cos-A LE  *  Kd  (8“10) 

where  Kd  is  a  partial-span  factor,  shown  plotted  in  Figure  69. 

The  method  for  generating  section  values  if  given  in  Section  8.4. 


Deflection  of  a  flap  produces  an  increase  in  lift  at 
zero  angle  of  attack  which  in  turn  produces  an  induced  drag 
given  by 


<*cLo)2 

CdI  =  KaKf  7TAR 


(8-11) 


where  Ka  and  Kf,  shown  plotted  in  Figure  70,  are  factors  which 
account  for  the  non-elliptical  span  loading  of  partial-span 
flaps . 


The  deflection  cf  a  flap  increases  the  camber  of  the 
airfoil.  In  Reference  37,  thin-airfoil  camber  theory  is  used 
to  relate  the  displacement  of  a  polar  with  flaps  to  the  lift 
increment  of  the  flap  at  zero  angle  of  attack  by  the  equation 

ZscL  =  (^cL)clean  +  CAcL)Due  to 

Wing  Flap 

where 

(ACOr..  - 


liiEL 


Flap  ^  1+1. 16^~~— j  ( . 5-Cf / C) 


(8-12) 


8 . 3  Moment  of  High-Lift  Devices 

The  pitching  moment  increment  caused  by  a  flap  on  a  swept 
wing  i^s  represented  by 

CM  -V  tanAc/2'AcJn  •  KSW  (8-13) 

oc  =0  oc  =  o  ° 

where  and  AC/  are  the  sectional  change  in  moment  and 

=»o  0 
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lift  at  zero  alpha  due  to  flap  deflection.  The  partial-span 
factors  Km  and  Ksw  are  shown  plotted  in  Figures  71  and  72,  and 
the  method  for  generating  sectiohal  data  is  given  in  Section  8.4. 


8 . 4  Aerodynamic  Characteristics  of  Two-Dimensional 
High-Lift  Devices 


The  lift  effectiveness  of  plain  trailing-edge  flaps  can 
be  estimated  from  thin-airfoil  theory.  The  rate  of  change  of 
lift  with  flap  deflection  at  a  constant  angle  of  attack  is 
given  by 


Ch 


-  2 


+  sin 


(8-14) 


where 

cos  ©£  ■  1  -  2(Cf/C*) 

and  C1  is  the  extended  chord. 

This  equation  is  plotted  in  Figure  73  as  a  function  of  flap- 
chord  ratio.  The  theory  considers  only  a  bent  flat  plate  and 
does  not  include  effects  of  thickness  or  large  deflection 
angles.  The  effects  are  accounted  for  in  References  38 
and  39  by  empirical  flap  efficiency  factors ;  as  shown  in 
Figures  74  through  77.  The  lift  of  a  plain  flap  may  now  be 
expressed  as 


A%  -  %  ■  c4  •  h 


(8-15) 


where 

is  the  plain- flap  efficiency  factor  from  Figure  74, 
depending  on  the  flap  deflection  angle, 
plus  the  included  angle  of  the  flap  trailing 
edge,  0TE. 

%  is  the  rate  of  change  of  lift  with  flap  deflection 
at  constant  angle  of  attack  from  Equation  8-14. 

$ f  is  the  flap  deflection  angle  in  radians. 


This  procedure  is  extended  to  slotted  flaps  with  Fowler  motion 
by  evaluating  C/c  at  a  flap- chord  ratio  based  on  the  extended 
chord.  For  double-  or  triple-slotted  Fowler  flap  segments  the 
lift  increment  is  obtained  by  summing  the  incremental  lift 
increments  for  each  flap  segment.  The  result  is 
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Figure  7  3  Theoretical  Lifting  Effectiveness  of  Trailing-Edge  Flaps 
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Figure  75  Turning  Efficiency  of  Single-Slotted  Flaps 
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(8-16) 


A%  =  Z1^SI±  '  C/S  t  '  ^fi 

i-1 

where 

i  is  a  subscript  that  indicates  the  1st,  2nd,  3rd  flap 
segment  of  the  slotted  flap 

I  is  the  number  of  slots  or  segments  in  the  flap  system 

^SJi  is  the  slotted-flap  efficiency  factor  from  Figures 
75,  76,  or  77  for  the  ith  flap  segment 

C/^  is  the  lift  effectiveness  for  the  ith  flap  segment 

Sft  is  the  flap  deflection  of  the  ith  flap  segment. 

The  method  of  summation  and  the  geometry  definition  required  to 
evaluate  Equation  8-16  is  shown  in  Figure  78. 

The  effects  of  leading-edge  high- lift  devices  on  the  wing 
lift  at  zero  angle  of  attack  is  estimated  from  thin-airfoil 


)  -  2(sin  eLE  -  eLE) 

LE 


(8-17) 


where 


cos  gLE  =  1  ~  2(Cle/C') 

Unlike  trailing- edge  flaps,  the  deflection  of  a  nose  flap  causes 

a  loss  in  lift  at  zero  angle  of  attack.  The  increment  in  lift, 

A  C/  is 
A  o 

4C,  -  <c,  )  .  £le  (8-18) 

o  lit/ 

where  S  hE  *-s  t*le  leading-edge  flap  deflection  angle  in  radians, 
positive  nose  down. 

The  two-dimensional  maximum  lift  increment,  A  due  to  a 

trailing-edge  plain-flap  deflection  is  given  in  Reference  38  as 

-  kt  ■  KS  ■  ^A-o  (4^)  (8-19) 
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and,  similarly,  for  single,  double-  or  triple-slotted  flaps  the 
equation  is 


Aca 


S.  K), 


AC4 

AC^ 


(8-20) 


where 


(4C'o). 


is  the  predicted  lift  increment  determined 
for  the  ith  flap  segment  from  Equation  8-18 


Kf  and  are  empirical  factors  developed  from  experi¬ 
mental  data  shown  plotted  in  Figures  79  and  30 

AC/ 

- . r - -  is  the  theoretical  relationship  between 

o  4C^max  and  AC/ ^  given  in  Reference  39  as 


(  max  \  _ 


9f 

Of+sinOf 


Jb\ j  sin|(X+9f)/sini(X-6f)| 


9f  tan  X/2 


(8-21) 


where  cosX  =  2(Xs/c')-l  and  cosGf  =  l-2(Cf/c').  This  equation 
relates  the  theoretical  maximum-lift  increment  to  the  chord  of 
the  flap  and  the  position  of  separation,  Xs/c',  on  the  airfoil. 

The  choice  of  the  separation  point  Xs/c’  to  determine  the  maximum- 
lift  ratio  from  Equation  8-21  depends  on  the  leading-edge 
Configuration.  For  clean  leading-edge  airfoils,  the  point  of 
flow  separation  is  assumed  at  the  leading  edge,  Xs/c’  =  0. 

For  airfoils  with  leading-edge  high- lift  devices,  the  point 
of  flow  separation  is  assumed  to  be  at  the  knee  of  the  leading- 
edge  device,  Xs/c'  *  C\&1  c' • 

The  two-dimensional  increment  in  maximum-lift  coefficient 
of  leading-edge  devices  is  predicted  in  Reference  38  as  follows 


‘  c^le)max-  7“* '  %  '  S] 

where, according  to  thin-airfoil  theory, 

C/  /-  )  =  2  sin  Orr 


(8-22) 


and  cose^g  =  (1-2Cee/c').  The  (C/^LE^max  presented  in  Figure  81. 
The  factors  ^7  max  and  are  empirical  factors,  introduced 

in  Reference  38*  to  correlate  A  C/max)LE  with  available  test 
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data  on  airfoils  with  trailing-edge  flaps  (Figures  82  and  83). 

The  rocxirouro-lif t  efficiency  factor,  Hmax)  depends  on  the  type  of 
leading-edge  device  and  on  the  ratio  of  leading-edge  radius  to 
maximum  airfoil  thickness;  is  an  efficiency  factor  that  accounts 
for  large  leading-edge-flap  deflections. 

In  the  case  of  the  two-dimensional  moment  increment,  the 
methodology  for  predicting  the  pitching  moment  is  developed  parallel 
to  the  methods  used  for  estimating  the  lift  increment,  which  extend 
thin-airfoil  theory  to  cover  multiple-slotted  flaps  with  extendable 
chords.  The  trail ing-edge-f lap  pitching-moment  increment  at  zero 
angle  of  attack  is  given  in  Reference  38  as 


^Cmoc=0  -  A  C/Q 

where 


(8-23) 


Ac/  is  the  predicted  lift  increment  for  either 

trailing-edge  or  leading-edge  devices 


l 


(CniA  is  the  theoretical  center-of -pressure  location 

j  from  thin-airfoil  theory  (Figure  84) 

kJ-j  is  an  empirical  factor  developed  from  experi¬ 

mental  data  (Figure  85) 


In  the  case  of  the  profile  drag  increment,  flap  drag  incre¬ 
ments  at  Ot  =0  for  plain  and  singie-slctted  flaps  are  obtained  from 
Figures  86  and  87.  These  figures  were  obtained  from  Section  6.1.7 
in  the  O.ATCOK.  For  double-  and  triple-slotted  flaps.  Figure  88  is 
used  to  obtain  the  =0  drag  increment. 


8 . 5  Ground  Effect 

During  takeoff  and  landing  when  the  clean  airplane  is  close 
to  the  ground,  the  ground  proximity  produces  an  increase  in  the 
lift-curve  slope,  ^  decrease  in  drag,  and  reduction  of  nose-up 
pitching  incraent.  However,  some  high-lift  configurations  may  show 
a  loss  in  lift-curve  slope  due  to  ground  effect. 
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Figure  84  Flap  Center  of  Pressure  Location  as  given 
by  Thin  Airfoil  Theory 
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Figure  85  Moment  Correlation  Factor  Versus  Flap  Chord 


I'igure  86  Two-Dimensional  Drag  Increment  Due  to  Plain 
Flaps 
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A  theoretical  approach  for  estimating  ground  effects  would 
employ  an  image*- vortex  system  to  represent  the  ground  plane. 

The  DATCOM  uses  a  semi-empirical  method  to  estimate  the  increase 
in  induced  upwash  at  the  wing  due  to  the  reflected  trailing 
vortices  and  the  change  in  lift  at  the  wing  due  to  the  reflected 
flap  ground  effects.  The  change  in  wing-body  lift  moment  and 
drag  due  to  ground  effect  at  a  constant  angle  of  attack  is  given 
by 


A  cl  “  +  (acl) 


(8-24) 


dCM  -  WCL)HTiHl/i 


.  acL2 

ACj>  "  ‘  7TAR  "  (%ing 

Alone 


rTCL 

57.3 


(8-25) 


(8-26) 


where  C r ,  r  and  T  are  ground- effect  parameters,  shown  plotted  in 
Figures  39,  90,  and  91-  The  term  Cjw  Alone  is  tbe  mini(nu,n 
drag  contribution  of  the  wing  plus  the  drag  due  to  lift  of  the 


wing  in  free  air.  The  increment  in  lift  on  the  wing  due  to 


ground  effect  is  calculated  from 


(t)  clm+(U)2«^0 


Flap 


(8-27) 


where  is  the  wing-alone  lift  at  angle  of  attack,  and  CT  , 

(r.L/L0),  and  A  (ACL)Flap  are  shown  plotted  in  Figures  92,  93, 

and  <94.  The  increment:  in  lift  on  the  horizontal  tail  due  to 
ground  effect  is  calculated  from 


(^l>ht  ~  Ut)  CLott  /  C1  '  He) 


(8-28) 


where is  the  change  in  downwash  at  the  tail  given  by 


A  £  =  6 


bef f  +  (HH  ~  h)2 
beff  +  <hh  +  H)2 


(8-29) 


where  £  is  the  downwash  at  the  tail  in  free  air,  H  is  the  height 
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b/2 


89  Prandtl’s  Interference  Coefficient 
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Figure  91  Parameter  Accounting  for  Variation  in 

Longitudinal  Velocity  with  Ground  Height 
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Figure  93  Parameter  Accounting  for  Ground  Effect 
Lift  Due  to  Bound  Vortices 
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of  the  'c/4  of  the  wing  above  the  ground,  Hjj  is  the  height  of  the 
5/4  of  the  horizontal  tail  above  the  ground,  and  beff  is  the 
effective  wing  span,  defined  as 


^ef  f 


CLm°  +  (AC^)Flap 

.  (AC^Flap 


(8-30) 


where  ClWo  is  the  wing  lift  in  free  air  without  flaps, 
&Lo>  p.-^ap  is  the  change  in  lift  due  to  the  flap,  and 


The  ratio  (b^/fo)  is  given  in  Figure  95  as  a  1  unction  of  taper  ratio, 
and  (bf/b^)  is  given  in  Figure  96  as  a  function  of  flap  span. 

The  increments  in  lift,  moment. and  drag  due  to  ground  effect 
are  calculated  at  each  angle  of  attack.  These  effects  are  then 
added  to  the  free-air  calculations;  so  that  a  trimmed  condition  in 
ground  effect  can  be  calcvflated. 
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9.  DATA  COMPARISONS 


This  section  presents  comparisons  between  predictions  using 
the  Large  Aircraft  program  and  test  data  from  various  sources. 

9 . I  Systematic  Wing  Study 

Comparison  of  the  predicted  lift  curves  and  drag  polars 
with  the  systematic  series  of  wing  reported  in  Reference  12 
are  presented  in  Figures  97  and  98.  Reference  12  reported  on 
the  tests  of  11  wing-body  models  covering  systematic  variations 
of  sweep,  thickness- to -chord  ratio,  position  of  maximum  thick¬ 
ness,  camber,  and  aspect  ratio.  The  study  was  conducted  over 
a. Mach  range  from  0.23  to  0.94.  The  predictions  were  made  using 
the  geometry  presented  in  Reference  12.  The  wetted  areas  of 
the  wing  and  body  were  computed  internally  by  the  Large  Air¬ 
craft  program  and  the  friction  drag  solution  was  for  fully  tur¬ 
bulent  flow. 

The  predicted  drag  polar  shapes  up  to  lift  coefficients 
slightly  greater  than  the  separation  lift  coefficient  comes 
close  to  matching  the  experimental  data  in  most  cases.  The 
predicted  minimum  drag  is  quite  close  to  experimental  levels 
for  most  cases  except  for  some  test;  data  below  0.3  Mach  number. 
Transition  grit  was  fixed  on  the  leading-edge  of  the  model  but 
its  possible  that  at  the  low  Reynolds  and  Mach  numbers  the  flow 
wasn't  fully  tripped  causing  the  test  minimum  drags  to  be  low. 
The  predicted  polars  in  the  transonic  region  tended  to  have  too 
much  drag  at  the  higher  lift  coefficients  probably  because  the 
drag  rise  contribution  was  increasing  too  fast. 

The  predicted  lift  versus  alpha  curves  are  in  good  agree¬ 
ment  with  experimental  data  at  the  low  lift  coefficients.  How¬ 
ever,  the  agreement  is  poor  in  the  high  lift  regions,  when  the 
predicted  maximum  lift  coefficients  don't  agree  with  test  data. 
At  the  low  speed,  low  Reynolds  number  condition  Cl^ax  was  under¬ 
predicted  especially  for  the  high  sweep  and  the  thin  wing  con¬ 
ditions.  The  thick  wing  lift  curve  slope  was  underpredicted 
at  the  highest  Mach  number  because  the  predicted  transonic 
thickness  correction  became  too  small . 


9 . 2  Cranked  Wing  Study 


Comparison  of  the  predicted  lift  curve  slope,  drag-due- 
to-lif t  factor,  minimum  drag  and  pitching  moment  slope  with  test 
data  for  three  of  the  wing  models  reported  in  Reference  40  are 
shown  in  Figures  99  through  102,  Each  wing  was  planer  and 
was  mounted  separately  on  a  cylindrical  body  of  revolution 
which  had  a  Sears-Haack  nose.  Model  1  was  a  59  degree  lead¬ 
ing-edge  triangular  planform  and  Model  2  and  3  had  two  straight 
line  leading  edge  segments  of  different  sweep  angles  (referred 
to  as  "cranked"  planforms) .  The  three  models  were  designed  to 
have  the  same  exposed  span  and  exposed  area. 


The  predictions  of  the  lift  curve,  drag  due-to-lift, mini¬ 
mum  drag  and  moment  curve  for  the  Model  1  wing  agreed  closely 
with  the  experimental  data  except  for  the  moment  prediction  at 
Mach  2.94  and  the  minimum  drag  at  the  transonic  conditions  and 
at  Mach  2.94.  The  predicted  lift  curve  slope  for  Models  2  and 
3  matched  the  experimental  data  at  the  low  subsonic  Mach  number, 
but  the  program  did  not  predict  the  large  increase  in  lift-curve 
slope  at  transonic  Mach  numbers  and  the  program  overpredicted 
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and  transonic  predictions  of  drag  due-to-lift,  for  Models  2  and 
3,  agrees  well  with  experimental  data,  but  the  supersonic  values 
are  underpredicted  because  the  lift  curve  slopes  were  overpre¬ 
dicted.  The  predicted  pitching  moment  curve  slope  was  in  fair 
agreement  with  the  test  data  for  Models  1  and  2  up  to  a  Mach 
number  of  1.6;  beyond  1.6  Mach  the  pitching  moment  was  over¬ 
predicted.  The  minimum  drag  prediction  for  Model  2  was  in 
fair  agreement  with  experimental  data  except  in  the  transonic 
region.  The  minimum  drag  predictions  for  Model  3  were  on  the 
average  14  percent  higher  than  the  test  results. 


9 . 3  C-141A  Flight  Test  Data 

Comparison  of  the  C-141A  flight  test  drag  polars  with  full- 
scale  corrected  wind-tunnel  polars  and  Large  Aircraft  program 
results  are  shown  in  Figures  103(a),  (b) ,  and  (c)  for  M=0.7, 
0.75,  and  0.775.  These  results  indicate  that  the  Large  Air¬ 
craft  program  predictions  are  in  good  agreement  with  the  flight 
data  points  at  the  higher  values  of  Cg  and  lower  values  of 
Mach  number.  Comparison  of  the  flight  test  versus  predicted 
drag-rise  characteristics  are  presented  in  Figure  104  as  drag 
variation  with  Mach  number  for  constant  lift  coefficients. 

These  results  indicate  that  the  reason  why  the  drag  polar  was 
underpredicted  at  the  higher  Mach  num’  ■►'s  and  lower  lift  co¬ 
efficients  was  due  primarily  to  overpreo^cting  the  drag  rise 
Mach  number,  Mcr. 


The  C-141A  flight  test  data  shown  in  Figures  103  and  104 
were  obtained  from  Reference  41.  The  flight  test  data  in 
Reference  42  are  corrected  to  an  equivalent  rigid-aircraft 
condition  trimmed  at  a  c.g.  location  of  0.25  mean  aerodynamic 
chord  at  a  Reynolds  number  of  55  x  10^/MAC.  The  corrected 
wind  tunnel  polars  shown  in  Figure  103  were  obtained  from 
Reference  42  where  fully-corrected  model  test  data  was  extra- 
poled  to  full-scale  Reynolds  number.  The  predictions  of  the 
Large  Aircraft  program  were  made  using  the  full  scale  C-141A 
geometry  presented  in  Reference  41.  The  program  predictions 
are  trimmed  at  a  0.25  MAC  c.g.  location  and  55  x  10^/MAC 
Reynolds  number  condition.  The  program  prediction  used  a 
7  count  (0.0007)  miscellaneous  drag  increment  in  the  drag 
buildup  to  represent  the  roughness  drag.  The  7  count  rough¬ 
ness  increment  was  the  same  as  used  in  Reference  42  to  scale 
model  data  to  full  scale. 


9 . 4  High  Lift  Configurations 

Figure  105  shows  the  lift,  drag,  and  pitching  moment 
predictions  for  a  wing-body  configuration  compare'  h  test 
data  for  a  clean  wing  and  a  partial  md  full-span  ...cngle- 
s lotted  flap  (Reference  43).  The  program  results  are  in  good 
agreement  with  the  clean  wing  lift  moment  and  drag.  The  incre¬ 
ment  in  lift  due  to  flap  deflection  is  underpredicted  at 
C'LMAX’  but  the  relative  predicted  differences  between  the 
full-span  and  partial-span  flaps  were  similar  to  the  test 
data.  The  increment  in  moment  due  to  flap  deflection  was 
underpredicted  for  both  the  full -span  and  partial-span  cases, 
while  the  drag  increment  for  flap  deflection  was  overpredicted 
only  for  the  full-span  condition. 
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105  Comparison  cf  Predictions  with  Test  Data,  Single-Slotted  Flap 


10.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  Large  Aircraft  Aerodynamic  Prediction  Program  provides  the 
aerodynamic 1st  with  a  quick-response  capability  for  evaluating  the 
aerodynamic  characteristics  of  arbitrary  large  aircraft  or  perform¬ 
ing  design  trade  studies.  The  Large  Aircraft  program  offers  a  sig¬ 
nificant  improvement  over  hand  calculations  based  on  handbook 
methods  because: 

1.  The  program  can  consider  more  complex  relationships 
between  geometric  and  aerodynamic  parameters  that 

may  be  neglected  or  not  feasible  in  hand  calculations. 

2.  Typical  hand  calculations  require  a  long  series  of 
intermediate  calculations  and  chart  lookups  that 
require  a  great  deal  of  time  to  perform  and  are 
prone  to  error. 

3.  The  program  always  performs  the  calculations  in  a 
consistent  manner,  whereas  hand-calculation  results 
may  vary  between  individuals  doing  the  same  cal¬ 
culations. 

Before  using  the  Large  Aircraft  program  to  analyze  an  arbitrary 
configuration,  the  user  should  familiarize  himself  with  the  methods 
and  operation  of  the  program.  If  a  similar-type  configuration  with 
test  data  is  available,  it  should  b*  evaluated  first  with  the  program 
so  as  to  establish  limits  for  credible  lift,  moment,  and  drag  pre¬ 
dictions  for  those  types  of  configurations. 

The  modular  construction  of  the  Large  Aircraft  program  will 
allow  subroutines  to  be  added  or  replaced  to  incorporate  new  and/ 
or  improved  aerodynamic  prediction  procedures  as  they  become  avail¬ 
able.  Future  improvements  to  the  Large  Aircraft  program  should 
extend  the  program  to  better  handle  fighter  aircraft  by  including  buffet 
predictions,  transonic  maneuvering  devices,  higher  Mach  numbers,  etc. 
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